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1.0  ABSTRACT 


A  12.75-inch  diameter  vertical  cylinder  was  subjected  to  large 
laboratory  waves  with  Keulegan-Carpenter  number  up  to  17  at  the 
local  force  transducer  (20  at  the  still  water  surface)  and  Reynolds 
numbers  up  to  2.7xl05  (3.2xl05  at  the  still  water  surface).  Local 
force,  local  pressures,  total  force,  and  the  elements  for  determin¬ 
ing  the  overturning  moments  were  measured.  The  purpose  was  to 
produce  a  high  quality  data  set  in  the  wave  regime  where  both  the 
velocity  and  the  acceleration -dependent  components  of  the  wave 
forces  are  important.  It  was  found  that  the  forces  were  more  accel¬ 
eration-dependent  than  originally  anticipated. 

The  local  forces  were  predicted  quite  well  witf  the  two-term 
Morison  equation  using  the  kinematics  from  the  measurements,  the 
linear,  and  the  stream  function  wave  theories.  The  best  predictions 
used  the  measured  kinematics. 


2.0  INTRODUCTION 


2.1  Background 

On  27  September  1982  a  contract  was  established  between  the 
Naval  Civil  Engineering  Laboratory  (NCEL),  Port  Hueneme,  California, 
and  Oregon  State  University  (OSU),  to  conduct  an  experiment  pertain¬ 
ing  to  wave  forces  on  a  12. /5-inch  diameter  vertical  cylinder.  The 
design  work  for  the  equipment  began  in  October  1982,  and  the  equip¬ 
ment  was  delivered  on  site  in  February  1983.  Ihe  laboratory  tests 
were  performed  in  the  period  from  12-26  April  1983  and  data  analyses 
were  completed  during  the  Fall  1983.  these  final  data  analyses  were 
made  possible  by  a  modification  to  the  original  NCEL  contract  with 
finding  from  the  Minerals  Management  Service  (MMS)  through  NCEL. 

The  Morison  equation  relating  wave  forces  on  vertical  cylinders 
is  given  in  the  following  equation.  Ihe  force  per  unit  of  cylinder 
length,  f,  is  expressed  by 

2 

F  »  Cd  ^-u(t)  |u(t)|  +  Cffl  u(t)  (2.1-1) 

where  Cd  and  Cm  are  the  empirically  determined  drag  and  inertia 
coefficients,  respectively  (Cm  *  1  +  Ca,  where  Ca  is  the  added  mass 
coefficient),  D  is  the  cylinder  diameter,  p  is  the  mass  density  of 
the  surrounding  water,  u(t)  is  the  horizontal  water  particle  veloci¬ 
ty,  and  u(t)  is  the  local  fluid  acceleration  au/at. 

Many  data  have  been  generated  from  experiments  by  many  investi¬ 
gators  to  determine  the  drag  and  inertia  coefficients  as  functions 
of  the  Keulegan-Carpenter  number,  K  (=  UmT/D,  where  Um  is  the  maxi¬ 
mum  water  particle  speed  and  T  is  the  wave  period),  and  the  Reynolds 


number,  R  (=  Umd/v,  where  v  is  the  kinematic  viscosity  of  the 
water).  However,  many  reportings  have  been  questioned  for  one 
reason  or  another  and  only  a  few  have  very  reliable  values  for  these 
force  transfer  coefficients.  One  reliable  source  is  from  planar 
oscillatory  flow  laboratory  experiments.  However,  they  are  limited 
because  the  flow  is  not  truly  that  of  waves.  Other  data  sources 
have  their  own  limitations.  For  example,  field  experiments  rely 
heavily  on  current  measurements  that  are  made  at  some  distance  from 
the  force  transducer  and  wave  theory  is  ultimately  used  to  extend 
the  measurement  to  the  force  transducer  location,  lo  our  knowledge, 
no  set  of  we  1 1 -documented,  thoroughly  instrumented  measurements  have 
been  made  in  a  truly  controlled  wave  condition  and  reported  on  in 
the  open  literature. 

Even  the  data  from  planar  oscillatory  flow  have  sone  scatter 
and  indicate  that  Eq.  1  may  be  improved  to  reproduce  well  the  forces 
acting  on  a  cylinder  within  certain  ranges  of  K.  Therefore, 
Sarpkaya  (1981)  suggested  the  use  of  a  so-called  four-term  Morison 
equation,  wherein  two  terms  are  added  to  Eq.  1  to  take  into  account 
energies  in  the  force  spectrum  that  occur  at  the  3rd  and  bth  harmon¬ 
ics  of  the  fundamental  wave  frequency.  The  coefficients  for  the 
added  two  terms  can  be  determined,  essentially,  from  the  drag  and 


inertia  coefficients  as  determined  from  the  two-term  equation. 

The  form  of  the  four-term  MUJS  equation  used  for  this  work  is 


n  n2  p  D  U~" 

F  =  cd2e-u(t)  |u(t)|  +Cm\-p  u(t)  +  — -yjz  £AC3+Bc3exp(Cc3a^ 

2A 

-1  /?  PDUm 

cos{3u>t-A  3+B^ 3exp ( 3a ) JJ  +  ~^[J?  ^Ac5+Bc5exp^Cc50t^  * 

•  cos  {Sot  -  A  ^  ^A<^  5+  5  ® C gf* ) ]}  (2.1*2) 

where  a  =  (K  -  12.5)2,  a  =  2tt/T,  K  =  l^T/D,  x  =  (C*  -  Cm)/KCd  (C*  is 


the  potential  flow  value--for  a  cylinder,  C*  =  2.9),  and 


Ac3  = 

0.01 

bc3  =  0.10 

Cc3  =  -0.08 

a*3  = 

-0.05 

B*3  "  -°*35 

C,3  *  -°-°4 

ac5  = 

0.0025 

Bj.g  =  0.053 

Cc5  •  -0.06 

A4>5  = 

0.25 

8,5  s  °'60 

8,5  *  -0-02 

Only 

two  terms 

were 

aJded  to  the 

Morison  equation  because 

the 

spectrum 

analyses 

of  the  residues 

from  the  U-tube  experiments 

(Sarpkaya, 

1981)  showed 

appreciable 

energy  from  only  these 

two 

higher  ordered  terms. 

2.2  Purposes  of  Test  Program 

The  major  purpose  of  this  program  was  to  provide  a  high  quality 
set  of  data  of  force  measurements  and  pressure  measurements  for  a 
vertical  smooth  cylinder  of  appreciable  size  (12.75-inch  diameter) 
in  real  waves  with  a  reasonable  range  of  Keulegan -Carpenter  and 
Reynolds  numbers  so  that  different  methods  of  analyses  can  be  exam¬ 
ined  (for  this  and  future  efforts).  Also  required  were  comparisons 
of  the  force  transfer  coefficients  from  different  methods  of  analy¬ 
sis  and  comparisons  of  forces,  using  the  differently  derived  coeffi¬ 
cients. 
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i 

2.3  Scope 

j  The  work  entailed  constructing  a  12.75-inch  diameter  test 

cylinder  equipment,  conducting  the  laboratory  model  tests,  recording 
and  analyzing  the  data,  and  production  of  the  reports.  The  test 
cylinder  equipment  was  of  all  new  construction.  Local  forces,  total 
forces,  and  the  elements  for  determining  the  overturning  moments 
were  measured.  In  addition,  eight  pressure  transducers  measured 
peripheral  time-dependent  pressures  around  the  center  of  the  local 

force  transducer.  Wave  water  particle  velocities  were  measured 

opposite  the  center  of  the  local  force  transducer,  in  line  with  the 
wave  axis  (crests).  Ihe  water  surface  profile  was  measured  in  line 
with  the  wave  axis  and  the  center  of  the  cylinder. 

A  minimum  of  26  runs  were  required  to  satisfy  the  requirements 
of  the  statement  of  work  (SOW).  However,  a  total  of  47  runs  were 
made  since  testing  time  did  permit  them.  Data  from  26  runs  were 
analyzed  as  required  in  the  SOW.  len  waves  in  a  sequence  were 

recorded  from  which  eight  were  later  extracted  for  the  analysis. 
The  eight  wave  periods,  beginning  at  any  time  in  the  sequence, 

allowed  for  seven  complete  peak-to-peak  waves  to  be  defined. 

The  force  transfer  coefficients,  Cd  and  Cm,  were  derived  with 
least  square  methods,  using  linear  and  stream  function  wave  theories 
and  the  measured  kinematics.  In  addition,  Cd  and  Cm  were  determined 
from  the  Wave  Project  II  d?"’,  where  Cm  is  constant  and  Cd  varies 
with  time  through  the  wave  cycle  (Hudspeth,  et.  al.  19/4),  and  from 
Sarpkaya's  U-tube  results  (Sarpkaya,  1981).  The  predicted  forces 
were  then  compared  against  the  measured  forces,  wherein  the  pre¬ 
dicted  forces  included  calculations  using  Eqs.  2.1-1  and  2. 


-6- 


In  summary,  the  analysis  included  12  combinations  of  time- 
dependent  force  determinations.  Using  the  kinematics  from  each  of 
the  linear  wave  theory,  the  stream  function  wave  theory,  and  the 
measured  kinematics,  the  time -dependent  forces  were  determined  for 
both  the  2-term  and  the  4-term  Mori  son  equations,  using  and  Cm 
determined  from  the  least  squares  method.  Wave  Project  II  data,  and 
the  U-tube  data. 


3.0  EXPERIMENT  SETUP  AND  PROCEDURES 


3.1  Wave  Research  Laboratory 

The  experimental  work  was  done  in  the  O.H.  Hinsdale  Wave 
Research  Laboratory  at  Oregon  State  University  (WRL).  The  overall 
length  of  the  wave  flume  is  340  ft.,  of  which  129  ft.  constitute  the 

length  in  which  tests  can  be  performed  well  beyond  the  evanescent 
modes  of  the  wave  board  and  in  front  of  the  toe  of  the  beach.  The 
location  of  the  test  cylinder  in  the  WRL  is  shown  in  Fig.  3-1. 
Good,  repeatable  waves  can  be  produced  ranging  from  a  high  frequency 
of  1.0  Hz  to  a  low  of  about  0.12  Hz.  For  these  tests,  the  wave  fre¬ 
quency  range  was  from  0.5  Hz  to  0.12  Hz.  Up  to  16  channels  of 
information  were  recorded  digitally  in  real  time  on  a  PDP  11  com¬ 
puter.  Other  channels  of  information  were  recorded  on  an  analogue 
tape  recorder  and  then  digitized  at  a  later  time. 

3.2  Experimental  and  Data  Acquisition  Equipment 

3.2.1  Test  Cylinder 

The  design  of  the  test  cylinder  required  that  it  be  stiff  for 
the  applied  loads  so  that  the  natural  frequency  of  the  lowest  mode 
would  be  relatively  high  with  respect  to  the  forcing  frequencies. 
From  impulse  response  tests,  where  the  cylinder  was  hit  with  a 
rubber  hammer  and  the  time  response  was  measured  by  strip  chart 
recorder,  the  local  force  transducer  demonstrated  a  natural  frequen¬ 
cy  of  8.2  Hz  when  submerged.  The  highest  wave  frequency  in  this 
test  was  0.5  Hz.  On  some  runs,  a  distinct  transverse  force  was 
experienced,  due  to  vortex  shedding,  that  had  frequencies  higher 


than  the  wave  frequency  by  up  to  approximately  three  times.  In 
addition,  other  turbulences  may  cause  transverse  forces.  The  ampli¬ 
tude  spectra  for  transverse  forces  in  Vol.  II  were  examined  for 
maximum  significant  frequencies.  An  approximate  value  of  1.5  Hz  was 
determined  as  a  conservatively  high  estimate,  even  though  the  ampli¬ 
tudes  therefrom  were  fairly  small.  Some  other  very  small  energies 
exist  at  5.0  Hz  and  at  8.0  Hz  for  the  center  force,  east  direction 
(probably  due  to  the  vibration  of  the  local  force  transducer  from 
near-impact  loads  in  steep  waves).  Thus,  the  largest  ratio  of 
importance  of  excitation  frequency  to  transducer  frequency  was  about 
1.5/8  =  0.19.  The  largest  ratio  of  fundamental  wave  frequency  to 
transducer  frequency  was  .5/8  =  .062.  Iherefore,  vibrations  of  the 
local  force  transducer  had  a  negligible  effect  on  the  measurements 
as  can  be  seen  in  any  reference  on  dynamics  for  the  response  curve 
of  a  one-degree-of-freedom  vibrator.  The  signals  were  also  filtered 
with  a  two-pole  Butterworth  filter  with  a  cutoff  frequency  of  8 
Hz.  However,  the  signals  were  recorded  filtered  and  unfiltered  in 
case  future  reference  is  to  be  made  to  the  unfiltered  signals. 

Another  design  criterion  was  to  size  the  test  cylinder  so  that 
the  component  parts  could  be  handled  relatively  easily  at  the  wave 
flume.  It  was  decided  to  provide  a  stiff  inner  core  structure  with 
a  removable  shell  surrounding  it. 

Figure  3-2  gives  the  general  concept  of  the  test  cylinder. 
Figures  3-3  and  4  are  photos  of  it  during  installation,  where  Fig. 
3-4  has  only  one-half  of  the  lower  skin  in  place.  Figure  3-5  shows 
the  strain  gauges  on  the  bottom  support;  Fig.  3-6  shows  the  bolt- 


-10- 


down  arrangement  at  the  bottom;  and  Fig.  3-7  shows  how  the  bottom 
skin  covered  the  supporting  core. 

The  upper  support  consisted  of  a  1.5-inch  round  stainless  steel 
rod,  30  inches  long,  inserted  into  roller  bearings  and  supported  at 
its  upper  end  with  a  U-joint  identical  to  the  one  at  the  bottom 
support,  as  shown  in  Fig.  3-2.  Thus,  the  structure  was  completely 
statically  determinate  and  no  longitudinal  forces  could  be  applied 
to  the  structure  from  temperature  variations  or  other  longitudinal 
stresses  from  construction  because  of  the  upper  roller  support. 

It  was  necessary  to  design  the  local  force  transducer  for  much 
smaller  forces  than  would  be  measured  at  the  top  and  bottom  trans¬ 
ducers  and  to  make  it  individually  statically  determinate  and  free 
from  the  small  deflections  of  the  central  structure.  Details  of  the 
local  force  transducer  (LFT)  are  shown  in  Figs.  3-8  and  9.  Photos 
of  the  strain  bar  supports  and  the  interior  of  the  LFT  are  given  in 
Figs  3-10  and  11. 

The  test  cylinder  was  supported  at  the  upper  end  by  an  8x12 
steel  tube  that  spanned  the  12  ft.  clear  distance  between  the  wave 
flume  walls,  to  which  it  was  securely  bolted.  The  wave  profiler  was 
mounted  in  this  beam  with  the  center  of  the  profiler  in  line  with 
the  center  of  the  test  cylinder,  perpendicular  to  the  wall.  It  was 
mounted  3  ft.  from  the  east  wall,  mid-way  between  the  wall  and  the 
cylinder  center  line.  A  Marsh-McBi rney  current  meter  was  mounted 
with  the  sensing  elements  directly  under  the  wave  profiler  at  a 
level  of  3.7  ft.  below  the  still  water  surface — the  same  level  as 
the  center  of  the  local  force  transducer.  A  second  Marsh-McBi  rney 


current  meter  was  mounted  directly  in  line  with  the  first  one,  but 
on  the  opposite  (west)  wall.  The  sensors  of  the  second  meter  were 
about  1.2  ft.  from  the  face  of  the  west  wall. 


The  test  cylinder  was  carefully  assembled  so  that  the  center 
line  was  vertical.  The  local  force  transducer  was  separated  only 
1/16-inch  from  the  rest  of  the  cylinder,  but  the  surface  was  slight¬ 
ly  out  of  alignment  with  the  smooth  cylinder  skins.  Figures  3-12 
and  13  show  the  best  and  worst  alignments.  Several  measurements 
were  made  of  the  misalignments  and  they  are  summarized  in  Table  3- 
l.  The  table  shows  the  LF I  was  slightly  larger  in  diameter  than  the 
adjoining  test  cylinder  by  0./  mm,  which  had  a  negligible  influence 
on  the  resulting  measurements. 

After  the  Druck  pressure  transducers  were  installed  and  the  LFT 
was  mounted,  the  surface  around  the  pressure  transducers  was  filled 
with  a  paint-like  substance  and  sanded  smooth  with  the  surrounding 
cylinder  surface. 

3.2.2  Force  Beam  Gauging 

Forces  were  measured  with  strain  gauge  force  dynamometers. 
Metal  foil  strain  gauges  were  used  of  120  ohms  each.  The  gauging 
stations  are  referred  to  as  the  top  gauge,  the  local  force  trans¬ 
ducer,  and  the  bottom  gauge.  Ihe  top  and  bottom  gauges  were  wired 
as  shown  in  Fig.  3-14. 

The  local  force  transducer  was  wired  as  shown  in  Fig.  3-15. 
Each  of  the  four  5/16-inch  square  support  beams  had  four  gauges, 
intending  to  have  two  complete  circuits  for  each  of  the  in-line  and 
transverse  directions.  However,  small  torsional  loads  from  the 


calibration  string  gave  inconsistent  results  so  the  two  in-line  and 
two  transverse  circuits  were  added  together.  This  action  nullified 
any  influence  on  the  readings  from  torsional  loads.  All  strain 
gauge  bridging  was  conditioned  with  Vishay  signal  conditioners. 

The  Vishay  strain  gauge  conditioner  contains  two  sections.  The 
first  is  a  precision  adjustable  voltage  regulator.  The  output  of 
the  regulator  is  used  as  the  excitation  for  the  strain  gauges.  For 
this  work  it  was  adjusted  for  10.00  volts.  The  exact  value  is  not 
critical  if  the  instrument  to  which  the  gauges  are  attached  is  cali¬ 
brated  at  that  voltage,  whatever  it  is.  The  second  section  includes 
several  stages  of  amplification  to  the  very  small  voltage  changes 
coming  from  the  strain  gauge  bridge.  Ihe  first  stage  is  a  differen¬ 
tial-type  amplifier  which  will  detect  and  amplify  only  voltage 
potentials  which  exist  between  the  two  return  leads  from  the 
bridge.  Any  potential  between  ground  and  one  lead  or  in  common  with 
both  will  not  be  amplified.  This  last  stage  then  provides  isolation 
from  many  noise  problems.  The  other  stage  of  the  anplifier  section 
provides  further  adjustable  gain  which  can  be  set  to  obtain  a  usable 
voltage  output  within  a  desired  range  from  the  conditioner  box. 

Auxiliary  circuits  in  the  conditioner  provide  bridge  balancing 
circuitry  to  nullify  any  static  offset  so  that  a  maximum  gain  may  be 
used  for  very  small  bridge  outputs.  Circuitry  for  different  types 
of  bridges  and  amplifier  adjustments  for  best  output  are  included. 
Oitputs  from  the  conditioners  were  input  to  the  Rockland  8.0  Hz  low 
pass  double  pole  filters. 


Selected  outputs  were  then  passed  into  summing-type  inverting 
operational  amplifiers  with  a  gain  of  very  nearly  -1.0.  The  exact 
value,  again,  does  not  matter  if  the  instrument  is  calibrated  using 
the  output  of  the  last  piece  of  electronic  equipment.  The  summing 
amplifier  has  a  very  large  bandpass  from  DC  to  several  thousand 
hertz  and  therefore  will  cause  no  measurable  phase  shift  at  our 
frequencies. 

Some  of  the  characteristics  of  the  general  operation  amplifier 
are  that  it  has  a  very  high  input  impedance  and  low  input  current 
(useful  to  prevent  loading  of  the  preceding  instrument).  Because  of 
this  very  high  input  impedance,  the  current  drawn  by  the  'op  amp'  is 
negligible  and  therefore  point  'A'  (see  Figure  3-16)  is  very  nearly 
at  ground  potential,  or  zero  volts.  Because  of  this,  the  current 
drawn  through  the  input  resistor  Rj  and  R2,  I.j,  is  equal  to  the 
current  through  Rf,  the  feedback  resistor.  Upon  solving  for  these 
relations  (without  including  every  step  here),  we  find  that  VQ/V^n  = 
-Rf/Ru.  Our  Rf/R^  is  set  to  be  1.0. 

In  the  more  specific  case  of  the  summing  op  amp,  the  current  1^ 
is  the  sum  of  the  currents  and  Ip2  and  so  we  arrive  at  the  rela¬ 
tion  (again  without  showing  all  the  intermediate  steps): 

-V0/Rf  «  V]/R1  +  V2/R2  (3.2.2-1) 

In  our  case  R^/Rf  *  1.0  and  R2/R^  =  1.0  and  so  we  have  as  a  final 
output  the  sum  of  the  inputs: 

-y0  =  vi  +  v2  (3. 2. 2-2) 
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This  being  the  last  instrument  we  have  our  final  output  voltage 
corrected  for  torsion  and  other  unwanted  bending  moments  and  with 


very  reduced  noise  reatty  to  be  digitized  by  the  computer. 

In  summary.  Fig.  3-17  shows  the  test  cylinder  with  only  one  of 
the  four  cylinder  skins  attached  to  the  core  beam.  Figure  3-18 
shows  three  of  the  four  skins  bolted  on.  Figure  3-19  shows  the 
completed  assembly  just  prior  to  testing.  The  west  wall  Marsh- 
McBirney  current  meter  is  shown  on  the  left  side  and  east  wall 
installation  is  shown  on  the  right  side  of  the  photograph. 

3.2.3  Data  Acquisition 

The  signals  were  routed  as  shown  in  Fig.  3-20.  A  maximum  of  16 
channels  of  information  can  be  digitized  simultaneously  with  the 
PDP-11  computer.  The  remaining  channels  were  recorded  on  analogue 
tape  for  digitizing  at  a  later  time.  The  analogue  tape  introduced 
additional  high  frequency  noise  so  the  signals  therefrom  were  again 
filtered  prior  to  digitizing. 

The  analogue  filter  is  a  two-pole  Butterworth  filter  manufac¬ 
tured  by  Rockwell.  rhe  digital  recording  was  done  with  a  PDP-11 
minicomputer.  The  analogue  tape  recorder  was  a  Bell  and  Howell 
model  CPR-4010.  rhe  digital  information  was  stored  on  an  il-track 
mag  tape  and  a  copy  of  all  data  recorded  were  provided  to  NCEL. 

In  order  to  help  synchronize  all  channels  in  real  time,  a  5- 
volt  step  function  signal  was  generated  by  the  computer  and  recorded 
on  one  channel.  This  5-volt  "switch"  was  activated  at  the  start  of 
digital  recording.  This  signal  was  also  recorded  on  the  analogue 
tape.  The  data  analysis  used  this  "switching"  technique  to  help 
phase  all  information. 


3.3  Test  Procedures 

3.3.1  Ca 1 i brations 

The  completed  cylinder  was  calibrated  in  four  force  directions 
corresponding  to  the  four  points  of  the  compass.  In  addition,  one 
horizontal  diagonal  calibration  was  made  to  illustrate  independence 
between  the  N-S  and  E-W  force  measurements.  Two  ranges  of  forces 
were  used  for  the  top  and  bottom  force  transducers  in  order  to  check 
their  linearity.  One  range  had  a  maximum  load  of  201.44  pounds  and 
the  other  was  16.56  pounds,  as  shown  in  Table  3-2.  Force  calibra¬ 
tions  were  made  at  one  or  more  of  the  five  elevations  as  shown  in 
Fig.  3-21.  One  measurement  was  made  with  the  forces  applied  simul¬ 
taneously  at  locations  2,  3,  and  4.  These  were  compared  with  force 
calibrations  at  position  3  alone.  The  purpose  of  this  was  to 
determine  if  distributed  loads,  which  create  bending  in  the  overall 
structure,  would  influence  the  force  measurements  at  position  3 
alone.  When  the  forces  were  applied  simultaneously  to  positions  2, 
3,  and  4,  it  is  designated  as  the  XXX  test.  In  that  case,  the  light 
calibration  was  applied  at  position  3  while  the  heavy  calibration 
loads  were  applied  to  positions  2  and  4  with  approximately  the  steps 
indicated  in  Table  3-2.  (This  loading  at  positions  2  and  4  could  be 
approximate  since  it  was  only  intended  to  provide  a  significant 
general  bending  of  the  core  beam  to  see  if  it  influenced  the  cali¬ 
brations  of  the  LFT.)  When  the  force  was  applied  to  position  3 
only,  it  is  designated  as  the  0X0  calibration.  It  was  found  that 
force  measurements  from  the  local  force  transducer  were  independent 
of  the  general  loading  on  the  cylinder.  Of  course,  during  tests  the 
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sensors  for  the  local  force  transducer  and  the  bottom  force  trans¬ 
ducer  were  submerged.  They  were  therefore  calibrated  in  the 
submerged  condition. 

The  XXX  calibration  arrangement  is  shown  in  Fig.  3-22.  Figure 
3-23  shows  a  calibration  in  progress  at  position  5  (submerged). 
When  calibrating  at  position  3  in  the  north  or  south  directions,  it 
was  necessary  to  provide  a  horizontal  truss  at  the  correct  level  to 
take  the  horizontal  and  vertical  thrust  from  the  calibration 
loads.  It  is  shown  in  Fig.  3-23. 

Linear  regression  analysis  was  used  to  obtain  the  calibration 
coefficients.  The  measured  calibration  data  values,  f,.  were  ap¬ 
proximated  by  a  linear  relation  to  the  measured  voltage,  x^, 
according  to 

f .  =  m  xi  +  b;  i  *  1,2,...,  N  (3. 3. 1-1) 

where  the  slope,  m,  and  the  intercept,  b,  were  determined  by  a  best 
least  squares  analysis,  fhe  correlation  coefficient,  R  was  computed 
by 

[z  fixi  -  z  f^x./N]2 

R2  =  --  « — - n — - 2 - 2 -  (3. 3. 1-2) 

[Z  r  -  (Zf.r/NJU  X.  -  (zx.)  /N] 

i  i  '  i  1  1 

and  is  a  measure  of  the  goodness-of-fit  of  the  two  regression  coef¬ 
ficients  since  0.0  <  R^  <  1.0.  If  R  were  identically  equal  to  1.0, 
the  data  would  be  perfectly  linear. 

The  results  of  the  force  calibrations  are  shown  in  Tables  3-3 
through  3-10. 
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The  low  range  of  calibrations  are  consistently  steeper  than  the 
high  range  of  calibrations,  but  they  only  differ  by  less  than  3.3% 
for  the  instruments  used  in  the  analysis.  The  calibrations  were 
quite  linear  and  consistent.  In  summary,  the  force  calibrations 
were  consistent,  linear,  and  very  repeatable  between  the  beginning 
of  the  tests  and  the  end  of  the  tests  for  all  of  the  data  channels 
used  in  the  NCEL  data  analyses.  Table  3-4  shows  that  TE1  Ch  #3  was 
not  consistent  and  should  not  be  used  in  future  data  analyses  of  the 


transverse  force  and  moments. 

The  cylinder  was  also  calibrated  at  location  5  with  the 
cylinder  in  torsion.  That  is,  the  forces  were  applied  to  the  edge 
of  the  cylinder  as  shown  in  Fig.  3-24.  This  was  done  in  both  the 
clockwise  and  counterclockwise  directions.  It  was  found  that  tor¬ 
sional  loads  had  very  little  effect  on  the  output  from  the  force 
transducers.  Therefore,  if  any  torsional  loads  were  applied  to  the 
cylinder  (and  they  could  only  have  been  so  from  uneven  boundary 
layer  development  between  the  east  and  west  sides  of  the  cylinder) 
they  would  not  influence  the  force  measurements. 

One  calibration  result  for  a  pressure  meter  is  shown  in  Fig.  3- 
25.  All  of  the  amplified  individual  calibrations  are  available  upon 
request.  They  were  all  very  linear  except  for  P6  (the  240° 

azimuth).  Figure  3-26  is  typical  of  the  linear  calibrations  and 
Fig.  3-27  shows  P6.  The  calibration  was  made  when  the  wave  flume 
was  being  filled  or  emptied.  Some  pressure  meters  had  a  consider¬ 
ably  larger  calibration  constant  than  others,  but  the  calibration 
values  remained  constant  and  the  pressure  measurements  were  consis- 


tent  and  reasonable  for  the  measurements  in  waves.  The  pressure 
calibrations  are  summarized  in  Tables  3-11  and  12. 

The  east  current  meter  had  been  previously  calibrated  by  Dibble 
in  (1980)  at  20.77  fps/V,  and  after  this  testing  by  Kim  in  June  1983 
at  20.47  fps/V  (a  1.5%  difference).  A  value  of  20.6  fps/V  was  used 
in  the  NCEL  data  analysis.  This  meter  was  more  consistent  than  the 
West  meter  and  was  the  one  used  in  the  data  analysis. 

3.3.2  Test  Plan 

The  test  cylinder  was  designed  and  constructed  in  the  Fall  and 
early  Spring  of  1982-83.  Prior  to  assembly  in  the  WRL,  the  force 
dynamometers  were  bench  tested  to  ascertain  their  workability.  The 
equipment  was  then  completely  assembled.  Calibrations  of  force 
transducers  were  made.  The  primary  current  meter  (East  side)  had 
been  previously  calibrated.  rhe  pressure  transducers  were  cali¬ 
brated  as  the  wave  flume  was  filled  and  emptied.  The  cylinder  was 
then  subjected  to  31  periodic  waves  to  satisfy  all  requirements  of 
the  SOW,  plus  a  few  extra  waves.  Then  four  random  wave  sequences 
were  generated  according  to  the  SOW.  Some  time  remained  in  that 
work  day  so  10  more  periodic  wave  runs  were  made  at  three  wave 
periods — to  fall  in  between  the.  wave  periods  previously  tested. 
(However,  only  26  periodic  wave  runs  were  processed,  according  to 
the  SOW). 

The  flume  was  then  pumped  down  and  calibrations  were  made  for 
the  force  transducers  for  comparison  with  the  calibrations  prior  to 
testing.  The  log  of  runs.  Table  3-13,  is  self-explanatory,  com¬ 
pleting  the  description  of  the  test  plan. 


Some  photos  of  the  cylinder 


Figs.  3-28,  29  and  30. 


4.0  DATA  ANALYSIS 


Periodic  wave  data  were  collected  on  22  April  1983  and  the 
random  wave  data  on  25  April  1983.  A  total  of  21  channels  of  elec¬ 
tronic  data  signals  were  recorded  during  the  NCEL  tests.  Table  4-1 
identifies  these  21  data  signal  channels  and  their  method  of  record¬ 
ing  during  both  calibration  and  wave  measurements.  Only  16  channels 
of  analogue-to-digital  (A/D)  were  available  on  the  PD  P-1 1/ E10  used 
at  the  WRL.  These  16  channels  are  identified  in  the  column  titled 
"A/D  Run  Ch.  #"  in  Table  4-1.  Ihe  remaining  five  channels  of  data 
recorded  during  wave  tests  were  recorded  on  a  Bell  and  Howell  CPR- 
4010  analogue  tape  deck  and  are  identified  in  the  column  titled 
"Analogue  Tape  Ch.#"  in  Table  4-1. 

The  local  force  transducers  9  and  10  were  combined  to  A/D  Cali¬ 
bration  Channel  #9  as  a  single  signal  for  both  digital  (A/0  Run 
channels  #5  and  #7)  and  in  analogue  (analogue  channel  #1)  format 
(see  Fig.  3-20).  A  five-volt  step  signal  was  used  to  reference  the 
analogue  time  series  to  the  discrete  digital  time  sequence.  This 
five-volt  step  signal  was  sent  by  the  PDP-11/E10  to  the  Bell  and 
Howell  CPR-4010  analogue  tape  recorder  at  the  instant  digital 
recording  began.  The  local  wave  force  transducer  signal  was 
recorded  simultaneously  on  two  A/D  channels  in  order  to  determine  if 
there  was  any  significant  phase  shift  introduced  by  the  8  Hz  Rock¬ 
land  filters.  The  A/D  Run  channel  #5  was  filtered  by  the  Rockland 
filter  while  A/D  Run  channel  #7  was  not  (see  Table  4-1). 

The  data  recorded  on  analogue  tape  were  passed  through  the 
Rockland  two-pole  Butterworth  filters  before  recording.  These 
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analogue  data  signals  were  again  passed  through  the  Rockland  filters 
before  being  digitized  in  order  to  filter  out  undesirable  high  fre¬ 
quency  noise  introduced  by  the  analogue  tape  recorder.  The  water 
surface  profiler  was  not  filtered  by  the  Rockland  filter  because  of 
signal  distortions  which  are  introduced  into  the  record  due  to 
"drop-outs"  in  the  data.  These  distortions  are  characteristic  of 
all  sonic  sensors.  The  question  of  phase-locking  the  data  properly 
to  the  water  surface  then  becomes  important.  All  the  data  records 
were  filtered  to  produce  the  proper  signal  amplitudes  and  phase 
angles  with  respect  to  the  water  surface.  The  theory  behind  these 
filters  is  presented  below. 

Six  channels  of  data  were  also  recorded  for  immediate  observa¬ 
tions  by  a  stripchart  recorder.  A  sample  record  is  shown  in  Fig. 
4-1. 


4.1  Record  Format  and  Filtering 

Figure  4-2  demonstrates  schematically  the  data  recording  pro¬ 
cedure  for  the  horizontal  current  meter  measurement. 

A  recorded  measurement,  m{t),  can  be  expressed  as  a  linear  sum 
of  cosines  and  sines  by 


m(t) 


=  a 

mo 


+ 


N 


l 

n=l 


(a  cos2tt f  t  +  b  sin2irf  t) 
mn  n  mn  n 


(4.1-la) 


where  m(t)  is  the  raw  (unfiltered)  recorded  value  of  the  horizontal 
current  measurement,  u(t),  and  amo  =  <m>,  the  average  value  of  m(t) 
during  the  recording  interval.  The  Fourier  transform  of  m(t)  is 
given  by 


Hm(f)  =  /  m(t  )exp(-i2irft)dt 


(4.1-2) 


where  the  complex-valued  Fourier  coefficient,  H  (f),  is  defined  as 


H  (f)  =  R  (f)  +  il  (f) 
m'  m  m 


or  HJf)  =  iMOjexptieCf)} 


(4  .l-3a) 
(4.1-3b) 


where  the  amplitude  (modulus)  is  defined  as 
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H  (f)|  =  (R„  +  IJ 
m  1  m  m 


(4  ,l-4a) 


and  the  phase  by 


i  I  (f ) 
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Substituting  Eq .  4.1-3b  into  Eq.  4.1-lb  demonstrates  that 


a  =  R  (f),  b  =  I  (f). 
m  m'  "  m  m'  ' 


(4.1-5) 


For  example,  high  frequency  noise  can  be  filtered  from  a 

measured  signal,  if  desired,  by  setting  a  =  b  =0  for  f  >  f 
3  3  mn  mn  c 

where  fc  is  the  "cutoff  frequency."  This  type  of  digital  filtering 
will  be  termed  the  “FFr  filter."  High  frequency  FFT  filtering  was 
done  to  many  of  the  NCEL  data  records.  For  example,  the  high  fre- 


quency  noise  above  1.5  Hz  was  filtered  from  the  horizontal  current 
and  force  measurements  (see  Fig.  4-10).  That  particular  cutoff 
frequency  is  specified  on  graphic  plots  when  it  was  used.  It  was 
selected  arbitrarily  to  be  low  enough  to  sufficiently  smooth  the 
signals,  but  high  enough  so  that  most  influences  from  vortex 
shedding  would  not  be  filtered. 

When  an  operation  is  made  on  a  signal  from  a  filter,  the 
Fourier  transforms  of  the  signal  and  filter  are  multiplied  in  the 
frequency  domain  in  order  to  determine  the  Fourier  transform  of  the 
output  from  the  filter.  If  the  signal  is  passed  through  two 
filters,  as  shown  in  Fig.  4-2,  there  are  two  multiplications.  Thus, 
the  Fourier  transform  of  the  recording  can  be  found  from 

Hm(f )  -  |Hu(f)|exp{ieu(f)}|^(f)|exp{ieM(f)})HR(f)|exp{ieR(f)} 

=  |Hm(f)|exp{iem(f)}  (4.1-6) 

where  1^  is  the  transform  of  u,  with  phase  0U  (both  quantities  are 
now  sought);  Hj^  is  the  transform  of  the  meter,  which  introduces  the 
phase  shift,  e^;  and  HR  is  the  transform  of  the  Rockland  filter, 
which  introduces  the  phase  shift,  0R.  To  simplify  the  following 
discussion,  it  is  assumed  that  |HR(f)|  »  1.0  and  0R(f)  B  0.  Later 
the  actual  values  can  be  considered  and  multiplied  by  the  product  of 
the  other  two.  Substituting  tqs.  4.1-3a  and  4.1-3b  into  Eq.  4.1-6 
gives 
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The  FFT  filtered  horizontal  current  u (t )  is  determined  from  the 
inverse  transform  according  to 


u(t)  =  /  H  (f )exp(i 2irft)df 


(4.1-9) 


Marsh-McBi rney  Filter 

The  single-pole  filter  characteristics  of  the  Marsh-McBi rney 
current  meter  are 


Vf)  =  1  +  i  2ttt  f 


(4.1-10a) 


-  [1  +  (Zirfr)'  ]  •  Cl  -  1  arft] 


(4.1-lflb) 


where  t  is  the  time  constant  of  the  current  meter  and  f  is  the 
cyclic  frequency  of  the  signal  in  Hz. 

The  comp  lex- valued  FFT  filtered  signal,  l^(f),  may  be  computed 


HJexp(i9m)  =  lHulexp(i0U)ltWexp(i0M> 


(4.1-11) 


which  may  be  used  to  compute  the  amplitude  of  the  current  from 


and  the  phase  from 


(4.1-12a) 
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The  value  of  the  time  constant  t  is  selectable  on  the  Marsh- 
McBirney  current  meters.  One  experimentally  determined  value  for 
the  east  meter  was  t  -  0.159  sec.  (see  Fig.  3  and  Table  4-1  in 
Dibble,  1980).  From  Eq.  4.1-10b  the  modulus  of  the  Marsh-McBirney 
filter  is 

|Hm|  =  [l  +  (2irtf  )^]"1/2  (4.1-13a) 

and  eM  =  TAN^t^nxf}  •  (4.1-13b) 

Table  4-2  summarizes  the  amplitude  attenuation  and  phase  shift  in¬ 
troduced  by  the  Marsh-McBirney  single  pole  filter  in  the  NCEL  data. 

The  current  meters  used  have  time  constants,  t,  with  three  or 
four  selections.  For  the  NCEL  tests,  they  were  set  on  0.2.  Small 
differences  in  phase  can  make  appreciable  differences  in  Cd  or  Cm, 
depending  on  whether  the  wave  is  drag  or  inertia  dominant.  The  time 
constants  used  for  the  NCEL  data  reduction  were  t  *  0.159  sec.  for 
the  east  meter  and  x  =  0.20  for  the  west  meter. 

Since  the  amplitudes  in  the  frequency  domain  are  proportional 
to  the  amplitudes  in  the  time  domain,  Eq.  4-13a  also  expresses  the 
division  correction  factor  for  the  amplitudes  of  each  frequency 
component  of  the  water  velocity  signal.  Likewise,  the  phase  angle 
properly  defines  the  phase  of  each  signal.  A  problem  was  encoun¬ 
tered  with  the  horizontal  signal  from  the  west  current  meter.  The 
source  of  this  problem  has  not  been  determined.  However,  the  east 
current  meter  (the  same  one  calibrated  by  Dibble,  1980)  gave  excel¬ 
lent  results  and  is  the  one  that  was  used  for  the  NCEL  data 
analyses. 
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Rock  land  Fi Iter 

The  Rockland  filter  is  a  two-pole  Butterworth  filter, 
filter  function  is  given  by 


(if)  +201foif+f;  (if)^2e2foif+fo 

where  fQ  is  the  cutoff  frequency,  f  is  the  signal  frequency  in  Hz, 
0!  =  0.924,  and  03  =  0.700  for  the  flat  decay  response.  Introducing 
a  dimensionless  frequency  ratio,  e  =  f/fQ,  reduces  Eq.  4-14  to 


HR(e)  * 


Me) 


(1-e  )+i  1.848e  (1-e  )+i  1.4e 

(1-4.5872  c2>e4)  -  i  3.248  E(l-e2) 


(1-4.5872  eS:  )  +  10.55  E^(l-2e^+e  ) 


(4.1-15a) 


(4 .1—1 5b ) 


The  amplitude  division  factor,  |Hr|,  is 

IHJ  =  L(l+1.415  e2>e4)(l-.04  e2+c4)j"1/2 


(4.1-16a) 


and  the  phase  angle,  0R,  is 


eo  •  tan-1  {- 


3.248  e  (1-e  )  i 
(1-4. 5872  e  +e) 


(4.1 - 16b ) 


Table  4-3  summarizes  the  amplitude  attenuation  and  the  phase  shift 
introduced  by  the  Rockland  filters  in  the  NCEL  data. 


Digital  Formatting 

Exactly  ten  periods  of  data  were  recorded  for  each  of  the  21 
channels  of  data  measured  during  the  periodic  wave  tests.  These 
data  were  reduced  to  exactly  eight  waves  and  were  digitized  at  a 


sampling  interval  determined  by  the  deterministic  wave  period,  T, 
and  the  FFT  base  2^  =  2048  according  to 

At  =  8T/2048  =  T/256  (4.1-17) 

Each  single  periodic  wave  in  the  eight-wave  sequence  thus  contained 
256  data  values.  This  data  digitization  scheme  was  required  in 
order  to  insure  that  the  data  could  be  FFT  filtered  as  previously 
described  using  an  FFT  algorithm  of  base  2.  The  FFT  algorithm  was 
taken  from  Bloomfield  (1976).  Only  seven  of  the  eight  waves  were 
actually  used  in  the  final  data  analyses. 

The  sonic  profiler  used  to  record  the  instantaneous  water  ele¬ 
vation  occasionally  suffers  from  data  ‘■drop-outs"  or  "over- 
ranges."  These  data  are  identified  as  "bad-points"  and  are  removed 
numerically  by  linearly  interpolating  between  the  good  data 
values.  On  rare  occasions,  these  intervals  of  bad-points  are  ex¬ 
cessively  long  and  result  in  forming  a  fictitious  wave  (see  Wave  #4 
in  Run  #5  in  Fig.  4-5).  These  data  have  been  ignored  in  the  data 
analyses  because  there  were  only  two  "bad  point"  waves  in  the  total 
of  182  waves  recorded.  These  256  data  points  per  wave  were  reduced 
to  33  data  values  (32  time  increments)  for  each  individual  periodic 
wave  in  the  final  data  analyses  for  the  force  coefficients.  The 
total  number  of  discrete  data  values  digitized  for  the  eight  waves 
in  the  41  periodic  runs  that  are  identified  in  Table  3-13  was 


NPTS  -  (41)  runs  •  (21)  .  (8)  .  (256)  ^ 


*  1,763,328  pts 

These  data  were  FFT  filtered  to  produce  the  records  that  were  used 
in  the  force  coefficient  analysis.  The  total  number  of  discrete 
data  values  used  in  the  final  force  coefficient  analyses  was 

NPTS  *  (26)  runs  (21)  (7)  (33)  ^ 

=  126,126  pts. 

For  the  four  random  wave  tests  listed  in  Table  3-13,  a  total  of 
8192  data  values  were  recorded  during  each  test.  The  total  number 
of  data  values  recorded  during  these  random  tests  was 

NPTS  =  (4)  runs  (21)  (8192) 

=  688,128  pts. 

4.2  Data  Processing 

The  16  channels  of  digitally  recorded  data  and  the  5  channels 
of  analogue  recorded  data  were  all  copied  in  digital  format  as  inte¬ 
ger  values  (0-4095)  on  magnetic  tape  by  the  PDP  1 1/El 0  at  the  WRt. 
This  data  tape  was  then  read  by  the  CDC  CYBER  73  main  frame  conputer 
at  Oregon  State  University  for  subsequent  data  processing  by  a  suite 
of  three  FORTRAN-coded  algorithms.  Figure  4-3  illustrates  schemat¬ 
ically  the  flow  chart  for  the  data  processing. 
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4.2.1  PROGRAM  WAVES 


The  integer  data  obtained  from  the  PDP  11/E10  digitizing 
process  at  the  WRL  was  converted  into  physical  units  using  the  ap¬ 
propriate  calibration  constants  listed  in  Chapter  3.  Figure  4-4 
illustrates  the  calibration  constants  entered  into  the  program  for 

run  # 5.  Figure  4-5  illustrates  a  sample  wave  profile  for  run  #5  in 

integer  (uncalibrated)  values;  and  Fig.  4-6  illustrates  the  same 

wave  profile  record  following  conversion  from  integer  values  to 
physical  units.  Note  the  “bad  point"  values  in  the  individual  waves 
#3  and  #4.  Figure  4-6  is  also  inverted  compared  to  Fig.  4-5  because 
larger  integers  (i.e. ,  higher  voltages)  indicate  wave  trough  values, 
so  a  negative  calibration  constant  is  required  (see  Fig.  4-4,  Column 
#1:  "WAVE  PROFILE"  =  -.5333). 

Following  calibration,  the  "bad  points"  were  removed  numerical¬ 
ly  by  linear  interpolation.  On  some  rare  occasions,  the  time 

interval  of  "bad  point"  values  was  so  long  that  these  data  were 
interpreted  by  the  numerical  algorithm  as  an  "inverted  wave"  and 
were  not  removed  by  the  subroutine.  This  rare  event  is  illustrated 
in  Figs.  4-5  and  4-6  where  the  "bad  points"  in  waves  #3  and  #4  have 
not  been  removed;  but  have,  instead,  been  converted  to  waves.  This 
rare  event  only  occurred  in  the  NCtl  data  during  Runs  #5  and  #24. 
Figure  4-7  illustrates  the  "bad  point  wave"  in  Run  #24. 

Figures  4-8  to  4-10  show  the  more  typical  results  obtained  from 
the  "bad  points"  removed  algorithm  for  Run  #12.  All  of  the  wave 
profiles  recorded  during  the  NCEL  tests  may  be  seen  in  Vol.  II. 


Next,  all  data  channels  were  Fourier  analyzed  by  a  FFT  algo¬ 
rithm  with  base  2.  The  FFT  algorithm  was  written  by  Bloomfield 
(1976).  One-sided  amplitude  spectra  were  computed  and  visually 
inspected  for  high  frequency  phenomena.  The  force  amplitude  spectra 
were  displayed  graphically  out  to  10  Hz.  Figure  4-11  illustrates  a 
transverse  force  amplitude  spectra  from  the  local  force  transducer 
for  Run  #30.  In  order  to  obtain  better  visual  resolution  of  the 
spectral  amplitudes  at  low  frequency,  all  spectra  were  plotted  out 
to  only  4  Hz.  Figure  4-12  illustrates  a  sample  horizontal  current 
amplitude  spectrum  for  Run  #b  at  this  expanded  scale  out  to  only  4 
Hz.  A  digital,  low-pass  FFf  filter  was  then  applied  to  the  ampli¬ 
tude  spectra  in  the  frequency  domain  to  remove  high  frequency 
noise.  Inversion  of  the  FFT  algorithm  synthesized  the  filtered  time 
series  for  later  analyses.  Figure  4-13  illustrates  a  filtered  time 
series  of  the  horizontal  current  for  Run  #b  with  a  cutoff  frequency, 
fc  =  1.5  It.  Figures  4-14  and  4-lb  illustrate  the  effect  of  the 
low-pass  FFT  filter  on  the  bottom  transverse  force  for  Run  #5. 

Horizontal  water  particle  accelerations  required  for  the 
Morison  equation  were  computed  from  the  measured  horizontal  current 
by  differentiation.  Ihis  differentiation  was  computed  in  the  fre¬ 
quency  domain.  Figure  4-16  illustrates  a  typical  horizontal 
acceleration  amplitude  spectrum  computed  from  the  measured  data  for 
Run  #5.  Although  the  horizontal  current  amplitude  spectrum  had  been 
FFT  filtered  at  l.b  Hz,  Mg.  4-16  illustrates  the  problem  of  ampli¬ 
fication  of  high  frequency  spectral  amplitudes  that  is  common  in 
numerical  differentiation  algorithms.  In  order  to  avoid  this  high 


frequency  amplification,  only  the  fundamental  and  the  first  harmonic 
were  used  to  synthesize  the  horizontal  acceleration  time  series. 
Figure  4-17  illustrates  the  horizontal  acceleration  amplitude 
spectrum  used  to  synthesize  a  horizontal  acceleration  time  series 
from  the  measured  horizontal  current.  Note  that  only  the  amplitudes 
at  the  fundamental  and  the  first  harmonic  have  been  retained  and 
that  ALL  of  the  other  amplitudes  have  been  set  equal  to  zero. 
Figure  4-18  illustrates  the  horizontal  acceleration  time  series 
synthesized  from  the  amplitude  spectrum  shown  in  Fig.  4-17. 

4.2.2  PROGRAM  FDATA5/FD5A 

The  FFT  filtered  data  from  PROGRAM  WAVE8  were  read  into  PROGRAM 
FDATA5  or  FD5A  for  the  selection  of  seven  individual  periodic  waves 
from  the  continuous  eight  wave  sequence. 

PROGRAM  FD5A  automatically  selects  the  first  seven  waves  in  the 
sequence  beginning  with  the  first  crest  as  shown  in  Fig.  4-19.  The 
period  of  each  individual  wave  in  the  seven-wave  sequence  is  deter¬ 
mined  by  first  establishing  a  time  at  a  wave  crest,  e.g.  crest 
number;  secondly  by  adding  one  wave  period  to  (as  set  into  the  func¬ 
tion  generator  that  drives  the  wave  board)  that  time.  This  accounts 
for  the  small  variability  in  wave  periods  found  in  the  data  analyses 
tables  in  Section  5.  However,  the  small  error  involved  averages  to 
zero.  The  wave  heights  are  confuted  by  averaging  the  leading  and 
trailing  crest  elevations  and  then  subtracting  the  trough  value  in- 
between  the  crests. 

PROGRAM  FDATA5  allows  for  the  manual  selection  of  the  individ¬ 
ual  waves  shown  in  Fig.  4-20.  Cross-hair  cursors  on  the  TEKTRONIX 
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4010  screen  are  positioned  at  the  individual  wave  crests  and  troughs 
to  define  the  seven  individual  waves.  The  wave  heights  and  periods 
are  then  computed  the  same  as  in  FD5A.  Figure  4-21  illustrates  the 
manual  determination  of  the  individual  waves  by  PROGRAM  FDATA5  for 
Run  #5.  Only  Runs  #5  &  #24  were  processed  by  PROGRAM  FDATA5  in 
order  to  manually  override  the  automatic  selection  of  the  fictitious 
"bad  point  waves"  that  would  have  occurred  had  PROGRAM  FD5A  been 
used. 

Finally,  both  programs  reduce  the  number  of  data  values  in  an 
individual  wave  from  approximately  256  values  (=  2048/8)  to  33 
points  per  wave  (32  time  increments).  This  reduces  the  CPU  time 
required  to  compute  the  force  coefficients,  Cd  and  Cm,  by  Fourier 
analysis  without  loss  of  accuracy. 

4.2.3  PROGRAM  LABT4 

PROGRAM  LABT4  computes  force  coefficients,  Cd  and  Cm,  from  the 
two-term  MOJS  equation  by  least-squares  analysis.  Each  of  the  seven 
individual  waves  selected  either  automatically  from  PROGRAM  FD5A  or 
manually  from  PROGRAM  F0ATA5  is  analyzed  separately.  First,  an 
average  velocity  is  computed  from  the  measured  horizontal  current 
profile.  Then  the  drag,  Cd,  and  inertia,  Cm  coefficients  are  com¬ 
puted  by  least-squares  analysis  according  to  Sarpkaya  and  Isaacson 
(1981).  The  predicted  2-term  Morison  equation  force  is  given  by 

F  (9)  =  1/2  CdpD[u(e)  +  <u>]|u(e)  +  <u> |  +  l/4Cmpir02  u(e)  (4.2-1) 


where  Fp(9)  is  the  instantaneous  predicted  local  in-line  force  per 
unit  of  the  vertical  cylinder  and  <u>  is  the  average  velocity.  A 
mean-square  error  between  the  instantaneous  measured  force  per  unit 
length,  Fm (0 ) ,  and  the  fore?  predicted  by  the  2-term  Morison  equa¬ 
tion  in  Eq.  4.2-1  may  be  defined  by 

s2-4r  /  CFm(8)  -  F  (e)]2de  (4.2-2) 
0  ■' 


The  drag,  Cd,  and  inertia,  Cm,  coefficients  may  be  estimated  from 
4.2-2  in  a  best  least-squares  sense  according  to 


3c  2 

air 

m 


=  0 


(4.2-3) 


By  Cramer's  rule,  the  estimates  for  the  two  force  coefficients  are 
given  by 


C  = 
m 


B1A4 


B2A2 


1/4  pirD  0ET{Cm,  Cd} 


C,  = 


B2A1  -  B1A2 
1/2  PD  0£T{Cm,  Cd} 


(4.2-3a) 


(4.2-3b) 


A1  =  /  u  (9  )de 

o 
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k2  -  f  u(e)[u(9)  +  <u>]ju(e)  +  <u>jde 
o 

r2lf  4 

A^  =  /  [u(9)  +  <u>]  d9 

o 

2tt 

=  /  Fm(e  )u(0  )d(9) 


where 


(4.2-4) 


b2  =  /  Fm(e)  Cu (0  )  +  <u>](u(6)  +  <u > | de 

o 

OET{Cm,  C,}  =  AxA4  -  A2 

The  integrals  in  Eqs.  4.2-4  were  evaluated  numerically  using 
Simpson's  rule  (Crandall,  p.  161,  1956). 

The  horizontal  water  particle  kinematics  required  in  the  inte¬ 
grals  in  Eqs.  4.2-4  were  computed  by  three  different  means:  1)  by 
the  measured  kinematics  from  the  east  current  meter;  2)  by  the 
linear  wave  theory  plus  the  average  measured  velocity;  and  3)  by  the 
nonlinear  Dean  stream  function  wave  theory  plus  the  average  measured 
velocity.  In  addition  to  the  three  sets  of  water  particle  kine¬ 
matics  used  to  obtain  the  force  coefficients,  two  additional  sets  of 
existing  force  coefficients  were  employed  in  the  analysis:  1)  Wave 
Project  II  coefficients;  and  2)  U-tube  coefficients.  Both  the  root- 
mean-square  and  the  maximum  errors  between  the  measured  and 
predicted  forces  were  computed  for  each  of  the  7  waves  analyzed  for 
each  Run. 

The  Wave  Project  II  force  coefficients  were  computed  by 


C  =  1.33 
m 


Cd(8>  “ 


'  1.20  ;  R  <  2.0xl05 

•  10E  ;  2.0xl05  <  R  <  5. 5xl06 

0.55  ;  R  >  5. 5xl06 


(4.2-5a) 


(4.2-5b) 


E  =  3.59197  -  1. 02271 og10(R)  +  0.0673774[1  og1Q(R) ]2 
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where 


R  =  j  u  (Q  )  +  <u  >  |  D/v 


and  where  R  is  the  instantaneous  Reynolds  number  and  v  =  kinematic 
viscosity  of  the  fluid  (=  1.41x10“^  ft^/sec. ).  Since  Cd  varies  with 
R  (and  time,  or  9),  there  is  no  single  value  to  show  for  a 
particular  run.  The  maximum  Reynolds  number  possible  during  any  run 
is  3x10^.  Therefore,  according  to  Eq.  4 . 2- 5b,  the  drag  coefficient 
will  vary  as  1.03  <  Cd  <  1.20. 

The  U-tube  coefficients  published  by  Sarpkaya  and  Isaacson 
(1981)  do  not  extend  to  either  the  Reynolds  numbers  or  the  Beta 
values  measured  for  the  NCEL  data.  Since  U-tube  tests  have  never 
been  extended  to  the  Reynolds  numbers  or  to  the  Beta  values  obtained 
in  large  wave  tank  experiments,  it  is  questionable  that  the  U-tube 
values  can  be  reliably  extended  to  in-service  design  of  prototype 
structures.  Nevertheless,  an  effort  was  made  to  extrapolate 
graphically  the  data  reported  for  U-tube  tests.  The  results  ob¬ 
tained  by  these  graphical  extrapolation  methods  are  not  unique  and 
should  not  be  interpreted  as  recommendations  for  the  use  of  U-tube 
coefficients  for  computing  wave  forces.  They  were  done  for  interest 
only  in  accordance  with  the  SOW. 

The  U-tube  inertia  coefficients  were  estimated  from 


C 


m 


'2.00  ;  K  <  2.0 

•  0. 001336K3-  0.029881K2+  Q.17002K  +  1.71668  ;  2.0  <  K  <  11.0 
1.75  ;  K  >  11.0  . 

(4.2-6) 


The  U-tube  drag  coefficients  were  estimated  from 
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R  <  l.OxlO5: 

Cd  =  0.000476K3  -  0.017877K2  +  0.185404K  +  0.244354  (4.2-7) 

R  >  l.OxlO5: 

r0.5  ;  K  <  2.5 

Cd  =  «  0.00124K3  -  0.3636K2  +  0.280324K  -  0.000779  (4.2-8) 

.0.55  ;  K  >  8.0 

Measured  Kinematics.  Figure  4-22  illustrates  the  measured 
horizontal  water  particle  kinematics  and  the  numerically  computed 
horizontal  water  particle  acceleration  for  wave  #3  in  Run  #5.  The 
measured  horizontal  water  particle  velocity  includes  the  average 
current,  <u>.  The  units  of  the  vertical  axis  in  Fig.  4-22  are  ft. 
for  the  water  surface  elevation,  ft/sec  for  the  velocity  and  ft/sec2 
for  the  acceleration.  Figures  4-23  to  4-25  compare  both  the 
measured  forces  and  the  forces  predicted  by  the  2-term  M0JS  equation 
using  the  measured  horizontal  water  particle  kinematics  and  each  of 
the  three  sets  of  force  coefficients. 

Linear  Wave  Theory  Kinematics.  Figure  4-26  illustrates  the 
horizontal  water  particle  velocity  and  acceleration  computed  from 
linear  wave  theory  for  wave  #3  in  Run  #5.  Figures  4-27  to  4-29 
compare  both  the  measured  forces  and  the  forces  predicted  by  the  2- 
term  M0JS  equation  using  the  linear  wave  theory  kinematics  and  each 
of  the  three  sets  of  coefficients  for  wave  #3  in  Run  #5. 

Table  4-4  summarizes  the  three  sets  of  force  coefficients  (Cd 
and  CJ  that  were  used  individually  with  each  of  the  three  different 
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wave  field  kinematics.  This  resulted  in  nine  different  analyses  for 
each  of  the  seven  individual  waves  in  each  of  the  deterministic  runs 
listed  in  Table  3-1  that  were  analyzed  by  the  2-term  MOJS 
equation.  (In  addition,  the  4-term  MOJS  equation  was  used  with  the 
three  sets  of  kinematics,  yielding  a  total  of  12  different  analy¬ 
ses.) 

Stream  Function  Theory  Kinematics.  Figure  4-30  illustrates  the 
horizontal  water  particle  kinematics  computed  by  the  stream  function 
wave  theory  for  wave  #3  in  Run  #5.  Figures  4-31  to  4-33  compare 
with  the  measured  forces  and  the  forces  predicted  by  the  2-term  MOJS 
equation  using  the  stream  function  theory  kinematics  and  each  of  the 
three  sets  of  coefficients  for  wave  #3  in  Run  #5. 

Four-Term  MOJS  Equation.  The  4 -term  MOJS  equation  developed  by 
Sarpkaya  (1981)  introduced  two  additional  time -dependent  terms 
having  universal  coefficients  that  could  be  used  with  the  Fourier 
coefficients  obtained  from  the  2-term  MOJS  equation  (Sarpkaya,  1981, 
p.  59).  Since  these  universal  coefficients  were  developed  by 
Sarpkaya  (1981)  specifically  for  the  use  with  the  Fourier  coeffi¬ 
cients  obtained  from  the  2-term  MOJS  equation,  it  was  not  realistic 
to  use  either  the  Wave  Project  II  or  the  U-tube  force  coefficients 
in  the  4-term  MOJS  equation  since  they  were  not  obtained  by  Fourier 
analysis  from  the  NCEL  data.  Consequently,  only  the  least-squares 
coefficients  were  used  for  each  of  the  three  sets  of  wave  field 
kinematics  in  the  analysis  by  the  4-term  MOJS  equation.  Sarpkaya 
and  Isaacson  (1981)  have  shown  that  the  inertia  coefficients  are 
identical  between  the  Fourier  and  least-squares  regression  analyses 


and  that  the  drag  coefficients  differ  only  slightly  (Sarpkaya  and 
Isaacson,  1981,  p.  90).  Therefore,  the  least-squares  force  coeffi¬ 
cients  were  used  with  the  universal  constants  in  the  4-term  MOJS 
equation  in  the  NCEL  data  analyses. 

Only  the  least-squares  coefficients  were  used  to  compute  the 
predicted  forces  by  the  4-term  MOJS  equation  using  each  of  the  three 
sets  of  kinematics.  Because  some  of  the  inertia  coefficients  com¬ 
puted  by  Eq.  4.2-3a  were  greater  than  2.0  and  some  of  the  drag 
coefficients  computed  by  Eq.  4.2-3b  were  negative,  the  lambda  term 
in  Eq  2.1-2  sometimes  becomes  negative.  When  this  occurred  the  two 
additional  terms  in  the  4-term  MOJS  equation  were  not  used. 
However,  if  the  inertia  coefficient  was  greater  than  2.0  and  the 
drag  coefficient  was  less  than  zero,  then  lambda  was  positive  and 
all  four  terms  were  computed. 

Figures  4-34  to  4-36  compare  the  measured  and  predicted  forces 
confuted  by  the  4-term  MOJS  equation  for  each  of  the  three  sets  of 
kinematics  using  only  the  least-squares  coefficients  for  wave  #3  in 
Run  #6. 

4.3  Random  Waves 

Four  random  wave  runs  were  recorded  during  the  NCEL  tests.  Two 
of  these  four  random  wave  runs  were  analyzed  and  10  individual 
periodic  waves  were  selected  for  analyses  based  on  comparisons  with 
wave  heights  and  periods  in  Table  3-13.  A  suite  of  three  conputer 
programs  was  developed  to  analyze  the  10  individual  periodic  waves 
that  are  listed  in  Table  4-5. 
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4.3.1  PROGRAM  WAVER 


PROGRAM  WAVER  calibrated  the  integer  data  from  the  random  wave 
runs  similar  to  PROGRAM  WAVE8  for  the  deterministic  wave  runs.  A 
block  of  data  containing  512  values  was  considered  as  a  sequence  of 
periodic  waves.  The  fundamental  wave  period  was  then  determined 
from  T  =  0.06*512  =  30.72  sec.,  where  At  *  0.06  sec.  This  single 
block  of  data  was  then  FFT  filtered  to  remove  the  high  frequency 
noise.  The  time  series  for  the  horizontal  water  particle  accelera¬ 
tion  was  computed  numerically  as  before  for  the  deterministic  wave 
data  in  the  frequency  domain.  However,  there  seemed  to  be  a  little 
more  high  frequency  noise  in  the  velocity  spectra  for  the  random 
waves.  Thus,  the  high  frequency  noise  problem  in  the  acceleration 
spectra  were  exacerbated  due  to  the  transfer  function  approach  that 
was  taken.  In  order  to  obtain  relatively  smooth  and  useful 
accelerations,  the  acceleration  spectra  were  generally  filtered  at 
0.5  Hz  instead  of  1.5  Hz,  which  was  the  general  case  for  the 
velocity  spectra.  Ihere  was  one  exception  for  Run  32,  wave  #2, 
where  the  horizontal  acceleration  was  filtered  at  0.4  Hz. 

4.3.2  PROGRAM  FDATAR 

The  FFT  filtered  data  from  PROGRAM  WAVER  were  then  visually 
inspected  on  the  TEKTRONIX  4010  CRT  screen  one  block  at  a  time. 
From  the  hard  copy  output,  the  individual  waves  were  selected  on  the 
basis  of  wave  height  and  wave  period.  Figure  4-37  illustrates  a 
single  block  of  512  values  of  the  wave  profile  from  random  wave 
#32.  The  individual  wave  #6  was  selected  for  analysis  as  an  indi¬ 
vidual  periodic  wave.  The  number  of  data  values  in  the  individual 
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periodic  wave  selected  was  reduced  to  33  values  for  that  period  by 
FFT  methods  (see  Appendix  3).  The  individual  waves  were  thus 
defined  from  crest-to-crest  from  the  water  profile  record. 

4.3.3  PROGRAM  LABT4 

The  individual  random  wave  was  analyzed  as  a  single  periodic 
wave.  Figures  4-38  to  4-52  illustrate  the  results  from  the  individ¬ 
ual  wave  #6  and  from  block  #7  in  random  wave  #32.  Note  that  in  Fig. 
4-38  the  wave  profile  is  asymmetric  and  that  both  the  sine  and 
cosine  components  were  retained  from  the  FF I  coefficients.  Ihe 
resulting  wave  profile  is  not  a  periodic  cosine.  Figure  4-38  also 
shows  that  the  velocity  profile,  as  measured,  was  not  perfectly  in 
phase  with  the  water  profile  as  measured  (which  was  from  crest-to- 
crest).  It  should  be  remembered  that  the  wave  period  was  defined  in 
terms  of  the  times  from  crest-to-crest  of  the  water  surface 
profile.  Because  of  current  measurement  error  the  measured 
horizontal  velocities  should  not  necessarily  be  in  phase  with  the 
surface  profile,  which  is,  in  fact,  the  case  for  Mg.  4-38. 

It  should  also  be  kept  in  mind  that  the  return  flow  (u)  is 
superimposed  on  the  water  particle  kinematics  when  a  periodic  wave 
theory  is  imposed  on  an  individual  wave  from  the  random  wave 
record.  therefore,  because  of  the  nonlinear  drag  term  in  the 
Morison  equation,  the  predicted  force  will  not  yield  a  symmetrical 
plot,  which  is  the  case,  for  example,  in  Fig.  4-43. 

Three  plots  of  the  comparison  of  wave  kinematics  appear  in 
Figs.  4-53,  54,  and  55.  Similar  plots  can  be  constructed  if  desired 
for  all  the  runs  from  the  results  shown  in  vol.  II. 
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5.0  SUMMARY  OF  RESULTS 


Table  5-1  summarizes  the  average  force  coefficients  from  the  26 
periodic  runs.  Table  5-2  summarizes  the  force  coefficients  from  the 
seven  individual  waves  for  each  of  the  26  periodic  wave  runs  and 
then  the  random  waves.  Each  of  the  seven  individual  waves  were 

represented  by  33  data  points  from  the  original  256  data  points. 
The  two  "bad  points'*  waves  in  Runs  #5  and  #24  are  identified  by 

asterisks. 

Figures  5-1  to  5-3  illustrate  the  parametric  dependency  of  the 
average  inertia  force  coefficient,  Cm,  on  the  Keulegan-Carpenter 
number,  K,  for  the  26  periodic  waves.  Figures  5-4  to  5-6  illustrate 
the  parametric  dependency  of  the  average  drag  force  coefficient,  Cd, 
on  the  Keulegan-Carpenter  number.  Note  that  negative  drag  coeffi¬ 
cients  are  not  plotted.  Also  shown  in  Figs.  5-1  to  5-6  are  the 
maximum  &  values  obtained  from  published  U-tube  tests  (Sarpkaya  and 
Isaacson,  1981)  and  the  3  values  for  these  tests. 

Figures  5-7  to  5-9  illustrate  the  inertia  force  coefficients 

for  the  10  individual  waves  selected  from  the  two  random  wave 

runs.  Figures  5-10  to  5-12  illustrate  the  drag  force  coefficients 
for  the  10  individual  waves  selected  from  the  two  random  wave 

runs.  Again,  negative  values  of  were  not  plotted. 

Table  5-3  summarizes  the  RMS  error  and  the  absolute  value  of 
the  maximum  error  for  both  the  26  periodic  wave  runs  and  the  10 
individual  waves  from  the  two  random  wave  runs.  The  RMS  errors  were 
calculated  from  the  differences  between  the  predicted  force  and  the 
measured  force  (as  functions  of  time).  In  each  case  the  predicted 


force  was  based  on  the  choice  of  kinematics  and  coefficients  as  so 
indicated.  For  example,  under  the  stream  function  theory  and  WP 1 1 
columns,  the  stream  function  kinematics,  in  conjunction  with  the 
force  transfer  coefficients  as  obtained  from  the  WP  1 1  data,  were 
used. 

Table  5-4  summarizes  the  RMS  error  as  a  percentage  of  the  maxi¬ 
mum  force  for  the  26  periodic  wave  runs  and  the  10  individual  waves 
from  the  two  random  wave  runs.  Also  tabulated  is  the  percentage  of 
energy  in  the  two-harmonic  acceleration  spectrum  compared  to  the 
unfiltered  acceleration  spectrum.  This  measure  indicates  the  loss 
of  energy  in  the  acceleration  spectrum  due  to  filtering  the  non¬ 
harmonic  components  in  the  spectrum.  Figures  5-13  to  5-21  display 
the  RMS  error  as  a  percentage  of  the  maximum  force  for  the  range  of 
Keulegan-Carpenter  numbers  measured.  These  results  were  obtained 
from  the  2-term  MOJS  equation. 

Figures  5-22  to  5-24  display  the  RMS  errors  as  a  percentage  of 
the  maximum  force  for  the  range  of  Keulagan-Carpenter  numbers 
measured  for  the  4-term  MOJS  equation.  Table  5-5  summarizes  the  RMS 
errors  for  the  analyses  dealing  with  the  4-term  MOOS  equation. 
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6.0  CONCLUSIONS 


Data  of  high  quality  were  obtained  on  measured  forces  and  wave 
conditions  for  a  12.75-inch  diameter  vertical  cylinder.  Drag 
and  inertia  coefficients  were  obtained  therefrom  for  wave  con¬ 
ditions  with  a  Keulegan-Carpenter  number  up  to  16  at  the  local 
force  transducer  and  20  at  the  still  water  surface.  According 
to  the  open  literature,  this  environment  should  have  produced 
conditions  that  were  drag  dominant  (although  the  inertia 
effects  should  still  be  significant).  However,  the  results 
indicate  the  conditions  were  inertia  dominant  (and  the  drag 
effects  were  still  significant). 

The  assumption  that  drag-dominated  wave  forces  exist  at 
Keulegan -Carpenter  numbers  near  approximately  20  is  not  cor¬ 
rect.  This  assumption  arises  from  U-tube  results  which 
demonstrate  that  there  is  no  parametric  dependency  of  the  force 
coefficients  for  K  >  20.  In  contrast,  the  Dean  number,  E,  (see 
Appendix  11.2)  demonstrates  that  wave  data  do  not  become  inde¬ 
pendent  of  K  until  E  >  4.0,  and  E  s  .088  K,  from  which  K  = 
45.  As  a  consequence,  most  laboratory  data  determined  herein 
are  dependent  on  both  Cd  and  C^  For  the  data  with  K  =  6  we 
found  that  the  center  forces  were  completely  inertia  dominated. 

The  wave  forces  as  functions  of  time  were  calculated  using 
measured  kinematics  and  theoretical  kinematics  from  linear  and 
stream  function  theory,  using  least  square  coefficients.  Wave 
Project  II  coefficients,  and  U-tube  coefficients--al 1  for  the 


2-term  MOJS  equation;  and  for  the  4-term  MOJS  equation  using 
least  square  coefficients.  The  3  values  were  much  higher  than 
normally  seen  for  laboratory  testing.  In  fact,  they  were  only 
moderately  smaller  than  for  conditions  at  sea.  The  Cm  values 
for  measured  kinematics  were  slightly  larger  than  for  the 
linear  and  stream  function  kinematics.  The  values  for  the 
measured  kinematics  had  less  scatter  than  for  the  two  theoreti¬ 
cal  kinematic  models.  In  addition,  the  values  of  in  the 

range  of  K  where  they  were  relatively  reliable  (K  >  10)  were 
approximately  one-half  those  reported  by  Sarpkaya  (1981)  for 
oscillatory  flow  U-tube  experiments  with  3  values  up  to  8370. 

3.  Plots  of  the  RMS  errors  vs.  K  indicate  that  the  most  accurate 
analysis  methods  used  herein  were  the  measured  kinematics  and 
the  coefficients  derived  from  a  least  square  analysis. 

4.  Ihe  published  U-tube  coefficients  do  not  extend  to  high  enough 
Reynolds  numbers  or  beta  parameters  to  be  used  for  the  OSU 
laboratory  wave  analyses.  Since  the  Reynolds  numbers  and  beta 
parameters  in  the  NCEL-OSU  laboratory  wave  tests  are  still 
lower  than  in-service  ocean  designed  values,  U-tube  coeffi¬ 
cients,  should  be  viewed  with  caution  for  wave  force-  design. 

5.  The  2-term  MOJS  equation  gave  lower  RMS  and  maximum  errors  than 
did  the  4-term  MUJS  equation.  ihe  universal  coefficients  in 
the  4-term  MOJS  equation  were  not  valid  for  these  laboratory 
wave  data. 
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The  2-term  MOJS  equation  used  with  measured  kinematics  and 
measured  coefficients  gave  lower  RMS  and  maximum  errors  com¬ 
pared  to  either  linear  wave  or  stream  function  wave  theory 
kinematics  or  the  Wave  Project  II  or  U-tube  coefficients. 

Linear  wave  theory  and  stream  function  wave  theory  kinematics 
should  be  corrected  for  wave  tank  circulation  currents  in  wave 
force  analyses. 

Although  interesting  observations  have  been  made  from  the 
several  types  of  calculations  required  by  the  SOW,  many  other 
analyses  wait  to  be  made  on  the  data.  Some  of  these  are: 

a)  Phase  shift  vs.  K  plots  for  water  surface  vs.  total 
force. 

b)  Maximum  force  coefficient  vs.  K  plots. 

c)  In  conjunction  with  a  and  b,  new  Cd's  and  Cm's  can 
be  derived. 

d)  Total  forces  and  overturning  moments  can  be  exam¬ 
ined. 

e)  Transverse  forces  should  be  examined  regarding  maxi¬ 
mum  lift  coefficients. 

f)  More  analyses  should  be  directed  to  the  random  wave 
runs.  Much  more  extensive  transfer  function  and 
wave-by-wave  analyses  can  be  made  to  shed  light  on 
proper  methods  and  techniques  that  differentiate 
analyses  in  random  waves  from  those  for  periodic 


waves. 


g)  For  both  periodic  waves  and  random  waves  a  cross¬ 
correlation  analysis  may  produce  more  insight  to  the 
phase  shift  model. 

9.  The  phase  shift  and  maximum  force  coefficients  analysis 
referred  to  above  has  innate  advantages.  First,  both  the  phase 
shift  and  maximum  force  coefficient  are  more  stable,  displaying 
little  scatter  for  a  1 1  values  of  K.  Secondly,  the  relative 

importance  of  the  drag  or  inertia  term  is  immediately  evident 
in  the  phase  plot.  thirdly,  the  maximum  force  is  much  more 
accurately  predicted  than  in  any  of  the  methods  utilized  in 
this  work.  It  should  be  added  that  this  proposed  technique  has 
only  recently  revealed  itself.  In  fact  the  ideas  encompassing 
the  new  methods  were  stimulated  by  this  research.  they  should 
be  further  explored  because  of  the  promise  of  a  major  improve¬ 
ment  in  available  methods  of  determining  Cd  and  Cm  from  labora¬ 
tory  experimentation  and  projecting  the  results  to  prototype 
conditions. 
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Table  3-1  Measurements  of  LFT  surface  alignment.  OUT  from  the  cylinder 
surface  =  +,  IN  =  Bearing  means  0°  =  North  (downwave) 

180°  =  South  (upwave)  and  90°  =  East.  Measurements  in  mm. 


Bearing 

(degree) 

8ottom 

A1 ignment 

Top 

A1  ignment 

0 

0.0 

+  .8 

30 

+  .3 

+  1.1 

60 

+  .4 

+  .8 

90 

+  .5 

+  .4 

120 

+  .5 

+  .2 

150 

+  .3 

0.0 

180 

+  .7 

0.0 

210 

+  1.0 

+  1.0 

240 

+  .8 

+  .9 

270 

+  .5 

+  .3 

300 

-  .6 

-  .5 

330 

-  1.1 

-  .5 

Longitudinal  Spacing 
From  Cyl i nder 

Top  Bottom 

1.4  1.5 


1.5  1.9 


1.5  1.3 


1.5  1.0 


Average 

a 


1.5  1.4 
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Ibs/volt 

correlation  coefficient 


Table  3-4  NCEL  VERTICAL  CYLINDER  CALIBRATION  SUMMARY 


UNITS:  lbs/vol t 


Table  3-6  NCEL  VERTICAL  CYLINDER  CALIBRATION  SUMMARY 


UNITS:  lbs/volt 

(•)  *  correlation  coefficient 


VERTICAL  CYLINOER  CALI 


UNITS:  lbs/vol t 

(■)  -  correlation  coefficient 


UNITS:  lbs/volts 

(•)  *  correlation  coefficient 


Table  3-9  NCEL  VERTICAL  CYLINDER  CALIBRATION  SUMMARY 


UNIT:  lbs/volt 

(•)  =  correlation  coefficient 


Table  3-10  NCEL  VERTICAL  CYLINDER  CALIBRATION  SUMMARY 


POST 


Table  3-11  NCEL  VERTICAL  CYLINDER  CALIBRATION  SUMMARY 
for  the  pressure  transducers 


‘correlation  coefficient 


NDER 


UNITS:  psl/v 

(•)  *correIation  coefficient 
(.{‘•value  fro«  UP-PRE  used  In  average 


fable  3-13  Log  of  runs 


Wave 

Target 

Anal og 

Counter 

Rema  rk  s 

Run  No. 

Frequency  T 

Span 

H 

Start 

Stop 

(Hz) 

(s) 

(wavehoard) 

(ft) 

1 

.500 

2 

95 

1.0 

35 

51 

2 

.500 

2 

32 

0.4 

51 

65 

X  10  A/D  .  2 

3 

.400 

2.5 

110 

1.1 

66 

83 

4 

.400 

2.5 

205 

2.0 

83 

101 

X  10  -  2 

5 

.400 

2.5 

375 

3.7 

101 

119 

X  10 

6 

.270 

3.7 

175 

1.0 

119 

140 

X  10  - 

7 

.270 

3.7 

444 

2.5 

140 

163 

X  10 

8 

.270 

3.7 

780 

4.4 

163 

185 

(voltage  out  of 

9 

.217 

4.6 

147 

0.7 

185 

209 

X  10 

10 

.217 

4.6 

336 

1.6 

209 

233 

X  10 

11 

.217 

4.6 

827 

3.1 

233 

257 

X  10  -  a 

12 

.217 

4.6 

1000 

4.3 

257 

281 

XI  -  #2 

13 

.189 

5.3 

202 

0.9 

281 

308 

X  10 

14 

.189 

5.3 

338 

1.5 

308 

334 

X  10 

15 

.189 

5.3 

693 

2.6 

334 

359 

X  10 

16 

.189 

5.3 

933 

3.5 

359 

386 

X  1 

17 

.189 

5.3 

1000 

3.9 

386 

412 

X  1 

18 

.167 

6.0 

665 

1.9 

412 

445 

X  10  4/25/83 

19 

.167 

6.0 

800 

2.6 

445 

477 

X  10 

20 

.167 

6.0 

950 

3.3 

477 

510 

X  1 

21 

.238 

4.2 

1000 

5.2 

510 

537 

X  1 

22 

.189 

5.3 

1000 

3.9 

537 

568 

X  1 

23 

.217 

4.6 

1000 

4.3 

568 

597 

X  1 

24 

.270 

3.7 

780 

4.4 

597 

622 

X  1 

25 

.400 

2.5 

375 

3.7 

622 

645 

X  1 

26 

.167 

6.0 

950 

3.3 

645 

676 

X  1 

27 

.5 

2 

95 

1.0 

676 

698 

X  10  Copy  of  R1 

28 

.270 

3.7 

780 

4.4 

698 

723 

X  1  Copy  of  R8 

29 

.217 

4.6 

827 

3.1 

723 

751 

X  1  Copy  of  Rll 

30 

.167 

6.0 

665 

1.9 

751 

784 

X  1  Copy  of  R18 

31 

.167 

6.0 

800 

2.6 

784 

816 

X  1  Copy  of  R19 

32 

Random 

JONSWAP  T  - 

3.7x  H1/3  .2.9 

816 

990 

X  1 

33 

Random  JON SWAP  T  * 

3.7x  Hj/3  «  2.9 

990 

1160 

X  1 

34 

Random  JONSWAP  T  - 

3*7x  Hl/3  *  2.9 

1160 

1330 

Oi  f  ferent 

35 

Random 

JONSWAP  T  » 

3.7x  H1/3  >2,9 

1330 

1498 

X  1 

36 

.333 

3.0 

140 

1.0 

1498 

1519 

X  1 

37 

.333 

3.0 

320 

2.5 

1519 

1539 

X  1 

38 

.333 

3.0 

850 

5.0 

1539 

1564 

X  1 

39 

.250 

4.0 

180 

1.0 

1564 

1587 

X  1 

40 

.250 

4.0 

450 

2.5 

1587 

1610 

X  1 

41 

.250 

4.0 

900 

5.0 

1610 

1633 

X  1 

42 

.200 

5.0 

200 

1.0 

1633 

1660 

X  1 

43 

.200 

5.0 

700 

2.5 

1660 

1684 

X  1 

44 

.200 

5.0 

1000 

5.0 

1684 

1710 

X  1 

45 

.333 

3.0 

650 

5.0 

1710 

1731 

X  1 

NOTES:  Water  Temperature  »  50^  (4/22  g  4/23) 
Water  Depth  -  11.5  ft. 


Table  4-1  Signal  routing  and  recording 


NOTE:  NR  »  Not  Recorded,  RL  =  Recorded  Later  (digitally) 


Table  4-2  Summary  of  amplitude  attenuation  and  phase  shifts 
introduced  by  the  Marsh-McBi rney  current  meters 
(t  =  0.159  sec.) 


fyl  9M 


(Hz) 

(radians) 

(degrees) 

0.0 

1.000 

0.0 

0.0 

0.2 

0.981 

-0.197 

-11.3 

0.4 

0.927 

-0.380 

-21.8 

0.5* 

0.895 

-0.463 

-26.5 

1.0 

0.707 

-0.785 

-45.0 

1.5 

0.555 

-0.982 

-56.3 

♦denotes  highest  wave  frequency  tested  during  periodic  wave  runs. 


Table  4-3  Summary  of  amplitude  attenuation  and  phase 
shifts  introduced  by  the  Rockland  filters 
(f0  =  8  Hz) 


(Hz) 

(radians) 

(degrees) 

0.0 

0.0 

1.0 

0.0 

0.0 

0.2 

0.025 

0.9996 

-0.080 

-4.6 

0.4 

0.05 

0.998 

-0.162 

-9.3 

0.5 

0.06  2  5 

0.997 

-0.202 

-11.6 

1.0 

0.125 

0.989 

-0.407 

-23.3 

1.5 

0.1875 

0.976 

-0.610 

-35.0 

-68- 

*  •**  *•  ,*»  *•  *. 
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Table  4-4  Summary  of  combinations  of  force  coefficients  and  wave 
kinematics  for  2-term  MOJS  equation  analyses 


Force  Coefficients 

1. 

Least-Squares 

(C(j  ar,d  Cm) 

2. 

Wave  Project  II 

3. 

U-Tube 

Wave  Kinematics 

1. 

Measured 

2. 

Linear  Wave  Theory 

3. 

Stream  Function  Wave  Theory 

Table  4-5  Individual  periodic  waves  selected  from  random  wave  runs 

Individual 
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Table  5-2  Summary  of  force  coefficients  for  each  of  the  seven  Individual  waves  In  the  26  periodic  wave  runs  and 
the  10  individual  random  waves  {‘denotes  “bad  points"  wave  omitted  from  analyses)  (RMS  err  *lQOx%) 
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vi  is  the  average  value  of  the  measured  horizontal  velocity  during  a  particular  wave.  Note 
that  it  is  mostly  negative,  indicating  a  flow  from  the  beach  toward  the  wave  hoard. 


Table  5-3  Suaatry  of  Maxlau*  end  RMS  Errors  for  2C  periodic  waves  end  10  Individual  rondo*  waves 
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1  Wave  Project  11  (Htdspetb,  et  el.,  1974) 
i  (Sapphaya  and  Isaacson,  1981) 
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Fig.  3-2  Sectional  elevation  of  test  cylinder 
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Fig.  3-5  Bottom  (and  top) 
support  strain  gaging. 


Fig.  3-6  Bottom  support 
bolting  to  WRL  floor. 


Fig.  3-7  Skin  installation 
at  bottom  support. 
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contiguous  when  assembled 
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Fig.  3-9  PORTAL  GAGE  FOR  LOCAL  FORCE 
TRANSDUCER  (exploded  view) 
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Fig  3-10  Local  force  transducer  strain  bar  support.  Pressure  ports 
shown  prior  to  filling  and  sanding  the  iMnedlate  area 
around  then. 


Fig.  3-11  Interior  of  the  LFT  with  8  pressure  transducers  and  4-1" 
round  bars  which  provide  structural  continuity  for  the 
core  structure. 
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Fig.  3-13 


The  worst  alignment  for  the  LFT  at  the  top  joint.  Offset 
*  1.1mm,  although  It  looks  larger  because  of  the  sun 
reflection. 
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Fig.  3-15  Local  force  transducer  strain  gage  wiring. 
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Fig.  3-17  Nearly  conpleted 
installation  with  one  of  the 
four  skins  attached. 


Fig.  3-18  Three  of  the  four 
skins  bolted  on.  Final  painting 
and  smoothening  around  the 
pressure  parts  in  progress. 


west  current  meter 


Fig.  3-19  Final  arrangement  just  prior  to  testing. 
Water  level  is  at  the  test  level. 


east  current 
meter 
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Fig.  3-22  XXX  Calibration  to  the  east  in  air. 


Fig.  3-23  Calibrating  at  position  5 
with  the  bottom  support  submerged. 
Also  shown  is  the  horizontal  truss 
used  for  calibrating  at  position  3. 
The  completely  assembled  vertical 
cylinder  is  on  the  left. 
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1  ?-  ftt-tj 


Fig.  3-26 Typical  pressure  transducer 
calibration  curve 


Fig. 3-27  P6  (240°)  pressure  transducer 
calibration  curve.  (Of  the  8  meters, 
this  is  the  only  nonlinear  one.) 


Fig.  4-3  Flow  Chart  of  NCEL  Data  Processing 
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CDC  CYBER  73 


PROGRAM  WAVES 


1.  Converts  integer  values  (0-4095)  into  physical  units  j 

2.  Removes  "bad  points"  from  data 

3.  FFT  data  with  fundamental  period  =  8T 

4.  FFT  filte*'  data  with  1.5  Hz  low-pass  digital  filter 

5.  Computes  horizontal  water  particle  acceleration 
spectrum  using  only1 2 3 two  harmonics 


_ i _ 

PROGRAMS  FDATA5/FD5A 

1.  Selects  7  discrete  waves  from  8  recorded  wave  sequence 
(FD5A  =  automatic  selection) 


(FDATA5  =  manual  selection  using  TEK  cross-hair  cursors) 

2.  Reduces  data  from  256  values/wave  to  33  values/wave 

3.  Produces  a  file  of  7  individual  waves  per  run _ j 


PROGRAM  LABT4 


1.  Computes  average  velocity  from  measured  horizontal 
water  particle  kinematics 

2.  Computes  C.  &  C  by  linear  regression  using  2-term  MOJS 
and  1)  measured?  2)  linear  wave  theory  +  average  velocity; 
and  3)  stream  function  theory  +  average  velocity. 

3.  Computes  errors  between  measured  and  predicted  in-line 
local  forces  using  2-  &  4-term  MOJS  with  3-sets  of  force 
coefficients  and  3  sets  of  wave  kinematics 
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Fig.  4-4  Sample  calibration  constant  format  for  Run  #5 
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Fig. 4-5  Sample  wave  profile  in  integer  values  for  Run  #5  (note  the 
"bad  point"  values  in  waves  #3  and  #4) 


05*21.00 


Sample  wave  profile  in  physical  units  following  calibration 
for  Run  #5  (note  that  profile  is  inverted  from  Fig. 4-5) 


fig.  4-7  Sample  wave  profile  with  "bad  point  wave"  for  Run  #24 


B3/1O/03.  01.31.54 


TIME  (SEC) 

Fig.  4-8  Sample  wave  profile  in  integer  values  for  Run  #12  (note  the 
absence  of  any  "bad  point"  values) 
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Fig.  4-9  Sample  wave  profile  in  physical  units  following  calibration  for 
Run  #12  (note  that  the  profile  is  inverted  from  Fig.  4-8) 


TIME  (SEC) 

Fig.  4-10  Typical  sample  wave  profile  without  any  "bad  points"  for  Run  #12 


83/10/tl.  08.04.50. 


Fig.  4-11  Sample  transverse  force  amplitude  spectrum  from  the  local 
force  transducer  for  Run  #30 
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Fig.  4-l£  Sample  amplitude  spectrum  for  the  horizontal  current  for  Run  #5 
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Fig.  4-13  Sample  FFT  filtered  time  series  for  the  horizonal  current 
for  Run  #5  (f  =  1.5  Hz) 


B3/10/04.  05-21.00 


Sample  time  series  of  the  bottom  transverse  force 
before  FFT  filtering  for  Run  #5 


TIME  (SEC) 

Fig. 4-15  Sample  time  series  of  the  bottom  transverse  force  after 
FFT  filtering  at  f  -  1.5  Hz  for  Run  #5 


Y*  Fig.  4-17  Sample  horizontal  acceleration  amplitude  spectrum  for 

Run  #5  (note  that  units  for  the  vert,  scale  should  read  ft/sec 


oenGRAM 


PROGRAM  FDATA5  for  Run  #5 


PROGRAM  FDATA5  STS^0”  0f  dividual 


04/01/16.  04.26.30.  T2SH37.005 


Fig.  4-24  Comparison  between  measured  (♦)  and  predicted  (■  )  local  in-line  force 
using  measured  kinematics  and  Wave  Project  II  coefficients  for  Wave  #3 
in  Run  #5 


B4/01/1B.  04.27.24.  T25H37.005 


between  measured  (♦)  and  predicted  (■)  local  in-line 
l  measured  kinematics  and  U-tub'e  coefficients  for  Wave  #3 


03/10/00*  00.47.17 


Sample  time  series  of  LWT  horizontal  water  particle  velocity 
and  acceleration  for  Wave  #3  in  Run  #5 


03/10/00.  OQ.40UB. 


-27  Comparison  between  measured  and  predicted  local  inline  force  using 
LWT  kinematics  and  measured  coefficients  for  Wave  #3  in  Run  #5 


4-28  Comparison  between  measured  and  predicted  local  inline  force  using 
LWT  kinematics  and  Wave  Project  II  coefficients  for  Wave  #3  in 
Run  #5 


03/10/00.  00.49^3 


Fig.  4-29  Comparison  between  measured  and  predicted  local  inline  force  using 
LWT  kinematics  and  U-tube  coefficients  for  Wave  #3  in  Run  #5 
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83/10/08.  00.52*27 


Fig.  4-32  Comparison  between  measured  and  predicted  local  inline  form 
using  stream  function  wave  theory  kinematics  and  Wave 
Project  II  coefficients  for  Wave  #3  in  Run  #5 


03/10/08.  00.53.11.  R25H37.C05 
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stream  function  wave  theory  kinematics  and  U-tube  coefficients 
for  Wave  #3  in  Run  #5 


Comparison  between  measured  and  predicted  forces  by  the  4-term 
MOJS  equation  using  measured  kinematics  for  Wave  #3  in  Run  #5 


03/10/08.  0050.12.  R25H37.00S 


Fig.  4-35  Comparison  between  measured  and  predicted  forces  by  the  4-term  MOJS 
equation  using  LWT  kinematics  for  Wave  #3  in  Run  #5 
(only  2  terms  used  since  C  >  2.0  and  C.  >  0.0) 


03/ 10/00.  00.KL56 


Fig.  4-36  Comparison  between  measured  and  predicted  forces  by  the  4-term  MOJS 
equation  using  stream  function  kinematics  for  Wave  #3  in  Run  #5 


83/10/21.  16.J/.-' 


83/10/25 


Fig.  4-38  Sample  tlmeseriesof  measured  horizontal  water  particle  velocity 
and  numerically  computed  acceleration  for  Wave  #6,  Block  #7t 
Random  Wave  #32 


Fig.  4-39  Comparison  between  measured  and  predicted  local  inline  force  using 
measured  kinematics  and  coefficients  for  Wave  #6,  Block  #7, 

Random  Wave  #32 


83/10/25.  8L34.31 


03/10/25.  0U37.32. 
UNCAR  WAVE  THEORY 


Comparison  between  measured  and  predicted  local  inline  force 
using  LWT  kinematics  and  Wave  Project  II  coefficients  for 
Wave  #6,.  Block  #7,  Random  Wave  #32 


03/10/25 


Fig.  4-46  Sample  time  series  of  stream  function  theory  horizontal  water 
particle  velocity  and  acceleration  for 
Wave  #6,  Block  #7,  Random  Wave  #32 


83/10/25.  01. 


Fig.  4-47  Comparison  between  measured  and  predicted  local  inline  force 
using  stream  function  wave  theory  kinematics  and  measured 
coefficients  for  Wave  #6,  Block  #7,  Random  WAve  #32 


03/10/25.  9LAU27 


Comparison  between  measured  and  predicted  local  inline  f 
stream  function  wave  theory  kinematics  and  Wave  Project 
coefficients  for  Wave  #6,  Block  #7,  Random  Wave  #32 


83/10/25.  0L.4Z.i-4 


Comparison  between  measured  and  predic 
stream  function  wave  theory  kinematics 
for  Wave  #6,  Block  #7,  Random  Wave  #32 


X  AXIS  IS  PHASE  IN  DEGREES 

4-54  Horizontal  velocity  comparisons  for  wave  #2  of  run  #29 


Horizontal  velocity  comparisons  for  wave  #2  of  run  #20 


PERIODIC  WAVE  -  STREAM  FUNCTION  THEORY 


Fig.  5-4  Average  drag  force  coefficient  for  the  26 
periodic  wave  runs  for  measured  kinematics 


Fig.  5-5  Average  drag  force  coefficient  for  the  26  periodic  wave 
runs  for  linear  wave  theory  kinematics 


PERIODIC  WAVES  -  STREAM  FUNCTION  THEORY 


Average  drag  force  coefficient  for  the  26  periodic 
wave  runs  for  stream  function  wave  theory  kinematics 


RANDOM  WAVES  -  MEASURED  KINEMATICS 


measured 


RANDOM  WAVES  -  LINEAR  WAVE  THEORY 


Fig.  5-0  Inertia  force  coefficients  for  the  10  individual  waves  selected 

from  the  two  random  wave  runs  using  linear  wave  theory  kinematics 


RANDOM  WAVES  -  STREAM  FUNCTION  THEORY 


Inertia  force  coefficients  for  the  10  individual  waves  selected  from 
the  two  random  wave  runs  usinq  stream  function  wave  theory  kinematics 


RANDOM  WAVES  -  MEASURED  KINEMATICS 


Drag  force  coefficients  for  the  10  individual  waves 
selected  from  the  two  random  wave  runs  using 
measured  kinematics  (negative  C.'s  not  plotted) 


Fig.  5-14  RMS  errors {%)  for  the  2-term  MOJS  equation 
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Fig.  5-16  RMS  errors  (%)  for  the  2-term  MOJS  equat 


Fig.  5-18  RMS  errors  (%)  for  the  2-term  MOJS  equation 


STREAM  FUNCTION  (LEAST  SQUARE) 


Fig.  5-19  RMS  errors  (%)  for  the  2-term  MOJS  equati 


Fig.  5-20  RMS  errors  (%)  for  the  2-term  MOJS  equation 
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Fig.  5-22  RMS  errors  {%)  for  the  •"♦-term  MOJS  equation 


LINEAR  WAVE  THEORY  <4-TERM  MOJS> 
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Fig.  5-23  RMS  errors  {%)  for  the  4-term  MOJS  equation 
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-24  RMS  errors  (%)  for  the  4-term  MOJS  equation 


11.0  APPENDIX 

11.1  Druck  pressure  transducer  specifications 

11.2  Parametric  dependency  of  force  coefficients  on  Dean  number 

11.3  Notes  from  Mr.  Standley  on  how  33  points  (32  time  increments) 
were  interpolated  for  each  random  wave  from  the  assigned  time 
increment  of  At  3  .06  sec. 


R  SPECIFICATIONS 


Druck  Limited 


3  Fir  Trw  Lana  Grobv 
Leicestershire  LE6  OFH  England 
Telephone  Leicester  (0533)  878551 
Telex  341 743 


Thousands  of  Druck  pressure  transducers  are  being 
used  throughout  the  world  in  applications  where  high 
accuracy  and  reliability  are  of  importance,  and  the 
PDCR 1 0  series  has  been  designed  more  specifically 
for  the  industrial  user  following  years  of  operation  of 
the  PDCR  22  and  PDCR  32  transducers  in  mainly 
aerospace  environments  requiring  maximum  per* 
formance  with  regard  to  weight,  reliability  and 
accuracy.  The  PDCR  1 0  is  an  extremely  rugged 
transducer  and  operational  ratings  such  as  shock, 
vibration  and  overload  are  extremely  high  as  is  usual 
with  microcircuit  devices. 

Pressure  ranges  extend  from  1  to  5000  psi,  gauge, 
sealed  gauge  or  differential.  The  standard  output  of 
1  OOmV  with  1 0  Volts  excitation  is  available  on  most 
pressure  ranges  whilst  still  retaining  high  iinearity  and 
hysteresis  accuracy  ( ±0-1%  B.S.L)  and  high  overload 
capability. 

A  silicon  strain  gauge  bridge  is  diffused  into  a  single 
crystal  silicon  diaphragm  using  microcircuit  techniques 
and  so  it  is  effectively  atomically  bonded,  and  this 
feature,  combined  with  the  fact  that  there  is  no 
measurable  hysteresis  in  such  diaohragm  material, 
makes  an  excellent  fundamental  sensing  element.  The 
gauge  factors  of  the  four  silicon  gauges  are 
approximately  50  times  those  of  conventional  wire 
strain  gauges  so  that  the  output  can  be  relatively  high 
whilst  still  retaining  a  high  overload. 


The  construction  of  the  transducer  is  completed  using 
microcircuit  alloying  and  bonding  techniques,  and 
electron  beam  welding.  The  ancillary  electronic 
package  to  give  good  temperature  compensation, 
rationalisation  or  amplification,  is  encapsulated.  The 
complete  transducer,  having  effectively  no  moving 
parts  and  being  both  mechanicailv  and  electronically 
extremely  rugged,  is  therefore  ideally  suited  to 
industrial  environments  of  extreme  vibration,  shock 
or  other  external  effects. 

Many  pressure  fittings  are  available  as  standard  and  the 
transducer  case  is  of  stainless  steel  and  titanium.  A 
polychloroprene  cable  is  sealed  into  the  body  with  a 
heat  shrink  boot.  The  integrity  of  the  cable  outlet  seals 
is  such  that  the  entire  transducer  could  be  immersed 
to  depths  of  at  least  300  metres. 

Different  versions  are  available  which  include  bore  hole 
transducers  suitable  for  continuous  immersion  in  salt 
or  high  mineral  water,  and  also  transducers  are  available 
with  a  2000  ohm  resistance  oridge  network  which 
makes  them  suitable  to  operate  with  two  wire  'mA' 
transminer  electronics. 

Also  current  outputs  can  be  sufficient  to  drive  meters 
and  galvanometers  of  osciilograoh  recorders  directly. 


Installation  details  for  PDCR 10 


Hexagonal  of  pressure 
connection  is  0-69'  A/F 
or  for  metric  1 7  mm.  A/F 


Standard  pressure  connections: 


B.S.P.  60*  inti,  cone 
t*  B.S.P.  60*  int».  cone 
■fa'  U.N.F.  as  MS.33656-4 


6din.M14x  1-5  mm 


Others  available  an  request 


Nominally  0-25  ins.  dia. 
4  core  P.C.P.  sheathed 


Scale:  actual  size 


Dimensions:  inches  (metric  mm.  in  blue) 


The  PDCR  10/D  BORE  HOLE  TRANSDUCER  has  the  same  body  dimensions  as  the  PDCR  lOabove,  except 
that  the  hexagon  and  pressure  fining  is  replaced  by  a  conical  front  end  as  illustrated  in  the  photograph  on  the 
previous  page. 

Electrical  Connection 

Red  -r  ve  supply.  White  -ve  supply.  Yellow  +  ve  output.  Blue  -ve  output. 

Ordering  information 

Please  specify  type  number  and  pressure  range  and  whether  this  is  gauge,  sealed  gauge  or  differential.  Also  state 
temperature  range  of  operation  and  survival,  and  degree  of  temperature  compensation  required. 

For  pressure  finings  refer  to  above  installation  drawing  but  where  non-standard  connections  are  required,  refer 
to  supplier.  The  cable  length  supplied  as  standard  is  1  metre,  and  longer  lengths  must  be  specified  for  which  an 
extra  charge  will  be  made.  Transducers  can  be  supplied  with  different  supply  voltages  and  different  outputs 
from  those  specified  and  in  all  cases,  reference  must  be  made  to  the  supplier. 

The  manufacturer  reserves  the  right  to  alter  specifications  without  notice.  &IQ  l  ^  ^  ^ 

Prices  also  subject  to  change  without  notice. 
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Druck  Limited 

3  Fir  Tree  Lane.  Groby,  Leicester  LE6  OFH, 
England 

Telephone:  Leicester  (0533)  878551 
Telex:  341 743 


Represented  by: 


SCANIVALVE,  INC 
p.  O.  BOX  20005 
SAN  DIEGO,  CALIF.  92120  U.S.A. 

TEL:  714  283-0010 
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Bore  Hole  and  Depth  Measurement 
Transducer  PDCR10/D 

A  version  of  the  PDCR 1 0  has  been  specifically 
designed  for  depth  measurement  in  bore  holes, 
reservoirs,  the  sea  and  many  other  applications.  It  has 
the  electrical  specifications  outlined  on  the  opposite 
page,  but  a  special  front  end  has  been  designed 
which  makes  it  ideal  for  dropping  down  small  bore 
pipes,  and  discourages  the  ingress  of  particles.  Extra 
integrity  has  been  built  into  the  cable  seal,  and 
combines  encapsulation,  0  ring  seals  and  a  heat  shrink 
boot.  The  pressure  ranges  are  vented  gauge  for  the 
lower  depths  but  sealed  gauge  versions  are  normally 
recommended  for  depths  above  about  1 50  feet. 

High  accuracy  versions  are  available  with  a  linearity 
and  hysteresis  combined  of  less  than  ±  0-06%  B.S.L 
and  temperature  compensation  from  — 2°C  to  +  30°C 
of  better  than  0-01  %/  °C  for  both  zero  and  span  shifts. 


Options 

For  many  years  we  have  specialised  in  producing 
pressure  transducers  based  on  the  use  of  an  integrated 
silicon  strain  gauge  diaphragm  and  have  solved  many 
different  pressure  measuring  and  recording  problems. 
Three  such  transducers  are  a  differential  version,  an 
amplifier  transducer  and  a  flush  mounting  transducer, 
and  these  are  briefly  mentioned  below. 

The  PDCR 1 0L  DIFFERENT1ALTRANSDUCER 
(Technical  note  4)  has  been  developed  to  fulfil  special 
applications  measuring  low  differential  pressure  at 
different  line  pressures  where  high  accuracy,  small  size 
and  fast  response  are  important.  These  transducers  have 
electron  beam  welded  pressure  connectors  at  each  end 
and  a  side  entry  cable.  The  current  models  are  suitable 
for  liquid  on  the  positive  side  and  dry,  non-corrosive 
gas  on  the  negative  side,  and  measurements  so  far 
have  been  limited  to  1 000  psi  static  pressures,  although 
other  types  are  being  developed.  The  zero  and 
sensitivity  change  for  increasing  line  pressure  is  often 
negligible  and  the  linearity  and  hysteresis  in  the 
negative  direction,  as  in  the  positive  direction,  is 
extremely  low. 

T ransducers  containing  I NTEG RAL  AM PU FI  ERS  are 
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also  being  manufactured  (Technical  note  5)  and  the 
more  standard  ones  require  a  power  supply  of  -i- 1 2.0, 
—1 2  Volts  and  give  output  voltages  of  up  to  5  Volts, 
with  a  low  output  impedance.  Again,  the  general 
performance  in  terms  of  accuracy  and  temperature,  are 
the  same  as  outlined  for  the  PDCR  1 0,  although  wider 
temperature  ranges  and  higher  accuracies  are 
available. 

The  FLUSH  MOUNTING  PRESSURETRANSUDCER, 
PDCR  50  (Technical  note  6)  has  a  0-5  ins.  diameter 
front  end  and  a  flange  which  enables  it  to  be  mounted 
into  tanks  and  other  containers  such  that  the  sensing 
diaphragm  is  flush  with  the  wall,  and  hence  gives  fast 
response  and  accurate  measurements  without 
disturbing  the  flow  pattern.  These  transducers  are 
already  proving  extremely  useful  in  marine  model 
testing,  general  flow  research  and  the  process 
industries. 

These  are  just  a  few  of  the  transducers  that  we  have 
manufactured  over  the  years  to  fulfil  specific 
requirements,  so  therefore,  if  you  have  pressure 
measuring  or  recording  problems,  it  is  suggested  that 
you  contact  us.  since  we  may  already  have  solved  your 
particular  problem. 


Pressure  Ranges*- 

± ‘1.2.5. 54050  1 00, 1 50. 250.  500. 1 000. 1 500. 

2500, 3500  5000.  !bsf/in  (p.s.i.). 

S.l.  units:  1 05N/m2  (bars)  0  to  6.895  kN/m1  through 
to  3.447  MN/mJ. 

Intermediate  ranges  are  available  on  request 

Operational  Mode 

Gauge,  sealed  gauge  or  differential. 

Overpressure 

400%  with  negligible  calibration  change  and  higher 
overload  capability  is  available  especially  in  lower 
ranges.  For  higher  ratings  please  specify. 

Pressure  Media 

Plus  side,  fluids  compatible  with  quartz,  titanium,  and 
silicon  dioxide. 

Input  and  Output  Impedance 
Nominally  1 000  ohms. 

For  optional  2000  ohms  refer  to  supplier. 

Excitation 

1 0  Volts  d.c.  regulated.  For  a.c.  excitation  please  specify. 

Output 
25mV  for  1  psi 
30mV  for  2.5  psi 
60mV  for  5  psi 
1 0OmV  for  1 0  psi  and  above 

Load  impedance  for  quoted  performance  should  be 
greater  than  100  Kohms.  For  lower  load  impedances 
please  specify. 

Resolution 

Infinite. 

Combined  Non-linearity  Hysteresis 
±0.1%  8. S.L 


Sensitivity  Setting 
±  3mV. 

More  accurate  settings  of  zero  and  span  are  available 
for  applications  where  interchangeability  is  important. 
Please  refer  to  supplier. 

Operational  Temperature  Range 
— 20°C  to  +  80°C  standard  but  this  temperature  range 
can  be  extended. 

Temperature  Effects 

Any  1 00°C  or  50°C  temperature  band  to  be  specified. 

Temperature  error  band  of  z0.5%  over  50°C 
or  temperature  error  band  of  r  1 .5%  over  1 Q0JC. 

Better  temperature  compensation  and  I  or  a  wider 
temperature  range  is  available  on  request.  Please  refer 
to  supplier. 

Natural  Frequency 

Typical  5000  p  Hz  where  p  is  full  range  pressure  in  psi. 

Acceleration  Sensitivity 

0.002%  FSO/g  for  5  psi  decreasing  to  0.0001 5% 

FSO/g  for  1 000  psi. 

Shock 

1 0OOg  for  1  ms.  in  each  of  three  mutually  perpendicular 
axes  will  not  affect  calibration. 

Weight 

Approx.  1 40  gms. 

Electrical  Connection 

Integral  sealed  cable  assembly:  standard  length  1  metre 
and  other  lengths  to  order. 

Mechanical  outline  and  ordering  information 
See  back  page. 


Zero  Offset 
±  3mV. 
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:e  coefficients  on  dean  number 


Ocean  Wave  Data.  Dean  and  Aagaard  (1969)  used  measured 
pressure  forces  from  a  prototype  offshore  oil  platform  and  a  least- 
squares  regression  analysis  with  theoretical  nonlinear  water 
particle  kinematics  computed  by  the  stream  function  theory  to  com¬ 
pute  Cd,  Cm.  They  found  that  the  drag  coefficient  varied  with 
Reynolds  number  and  that  the  inertia  coefficient  was  essentially  a 
constant.  The  Dean  and  Aagaard  regression  analysis  illustrates  the 
parametric  dependence  of  the  coefficients  and  explains  some  of  the 
scatter  observed  by  various  investigators. 

Dean  and  Aagaard  (1969)  computed  Cq,  Cj  from  a  least  squares 
fit  to  the  measured  forces  by  minimizing  the  following  mean  square 
error: 


£  *.<■>-%  <"»2 


Minimizing  this  mean  square  error  requires  that 


f§- •  0  .nd  f£- •  0  (2) 

m  d 

Equation  1  for  the  mean  square  error  may  be  expanded  into  a  general 
equation  for  a  quadratic  surface  of  second-degree  with  an  origin 
which  has  been  translated  and  rotated.  The  coefficients  of  the 


independent  variables  include  products  of  the  wave  field  kinematics 
and  measured  forces  and  determine  whether  or  not  the  data  are  well 
conditioned  or  ill  conditioned  for  determining  true  minima  (Dean, 


Expanding  the  mean  square  error  yields  the  following  equation 
for  a  quadratic  error  surface: 

e2  =  (“2“)  <u4>  +  (Cm  ^)2  <4j2>  +  A-J- ’)  (Cm  <u  | u  | u > 

-  2  (-P)  <Fm  UM>  -  2  (C.  ^)<V>  '  <Fi>  (3) 

in  which  the  temporal  averaging  operator  <•  >  is  defined 

i  !? 

by  <•  >  =  u  l  (*)n-  Th i s  equation  is  the  general  equation  for  an 
n=l 

ellipse  whose  origin  has  been  translated  and  rotated.  For  data 
which  are  simple  harmonic  oscillations,  the  coordinates  of  the 
origin  (xQ,  yQ)  are  given  by 

,  ,16,  <FmliM> 

x°  ( 3p )  4  (4a) 

^  rn 

J'o  =  C-%J  ((Tt)2  0  <FmS>  (4b) 

pir  m 

To  determine  the  suitability  of  the  data  to  resolve  the  regres¬ 
sion  coefficients,  Cd,  Cm,  changes  in  the  minimum  values  of  these 
coefficients  for  a  given  mean  square  error  are  minimized. 

r0r  Sn  =  -onstant,  it  may  be  shown  that  this  minimization  of 
the  changes  in  the  minimum  values  gives 


_/Z  -  (=f)  ">4"31/2 


and  for  C,  =  constant: 


The  ratio  of  the  axes  of  the  error  ellipse  (eccentricity)  is  given 
by 


E 


_  6Cm  _  2_  r<u4>a/2 
’  6Cd  '  *D 


For  simple  harmonic  oscillations,  this  reduces  to 

e  •  hj  wrt'1  f 

.  r,rtll/2  Cm  |Fd)max 


(6) 


(7a) 

(7b) 


The  error  surface  for  simple  harmonic  oscillations  in  dimensionless 
form  is 


2  2 

£  .  -  (c - )mi n  5C.  F*  6  C  n 

,r  x . ■  w  (aj  *  i/2  tsft 

(I-q)  max  m  D  d 


(8) 


Dean  (1976)  gives  the  following  criteria  for  estimating  the  suit¬ 
ability  of  the  data  for  determining  Cd,  Cm: 


E 

- 075 

0.25  -  4.0 
4.0 


Relatively  Vie  11 -Conditioned  to  Determine 

- - ^ - 

m 

C ,  and  C 
d  m 


From  Eq.  7a,  it  is  obvious  that  the  Dean  number,  E,  is  proportional 
to  the  Keulegan-Carpenter  number,  K.  The  constant  of  proportion¬ 
ality  for  simple  harmonic  flows  is 


E/K 


E  x  11.4 


(9) 


APPENDIX  11.3 


FDATAR 

Choose  starting  peak  of  desired  wave. 

DT  =  .06  sec 

NP  =  #  of  points  in  original  wave  profile  including  both  end  points 
at  spacing  of  DT 

PERIOD  =  (NP-1) *DT  period  of  desired  wave 


NFT  =  #  of  points  to  FF r  including  desired  wave  in  center 

OUTPUT  wave  must  have  NN  =  33  points 

Spacing  of  output  wave  (delta-t)  =  IT  =  PERIOD/ (NN- 1) 

N  =  #  of  spaces  in  total  FFT  window  at  delta-t  =  TT 

FFT  the  data  and  filter  at  2  Hz  since  original  data  were  already 
filtered  at  1.5  Hz 

Starting  at  time  -  LAG*DT  compute  NN  points  using  the  coefficients 

NFLTR-1 

WAVE  ( I )  =  XR  ( I )  +  2.*  I  XR(m)*cos(2n*(m-l)*T/(N*TT) ) 

m=2 

-  XI (m)  *si n (2ir*(m-l)*r/ (N*TT) 

where  T  =  (LAG*DT)  +  TT*(I-1). 
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By  direction 


Rockland  Fi Iter 


The  Rockland  filter  Is  a  two-pole  Butterworth  filter.  The 
filter  function  Is  given  by 


WO 


Of)^282f01f+fJ 


(4.1-14) 


where  f0  Is  the  cutoff  frequency,  f  Is  the  signal  frequency  In  Hz, 
*  0.924,  and  $2  *  0.700  for  the  flat  decay  response.  Introducing 
a  dimensionless  frequency  ratio,  e  a  f/fQ,  reAices  Eq.  4-14  to 


Me) 


(1-0+1  1.848c  (1-0+1  1.4e 


_  (1-4.5872  e^+e4)  -  1  3.248  e (1-e2) 

HD(e)  *  J - 5 — rt - 7~ - 9  ~r 

R  (1-4.5872  t+tf  +  10.55  eZ(l-2eZ+e) 


The  amplitude  division  factor,  jHR|,  is 


1^1  -  1(1+1.415  eZ*e4)(l-.04  e2+e4)j"1/2  (4.1-16a) 

and  the  phase  angle,  eR,  Is 

9„  ■  TAN'1  {-  3.248  ,  (4.1-16b) 

R  (1-4.5872  c  ) 

Table  4-3  summarizes  the  amplitude  attenuation  and  the  phase  shift 
Introduced  by  the  Rockland  filters  In  the  NCEL  data. 


(4.1-15a) 


(4 .1— 15b ) 


(4.1-16a) 


Digital  Formatting 

Exactly  ten  periods  of  data  were  recorded  for  each  of  the  21 
channels  of  data  measured  during  the  periodic  wave  tests.  These 
data  were  reduced  to  exactly  eight  waves  and  were  digitized  at  a 


sampling  Interval  determined  by  the  deterministic  wave  period,  T, 
and  the  FFT  base  2^  »  2048  according  to 

At  -  8T/2048  -  T/256  (4.1-17) 

Each  single  periodic  wave  In  the  eight-wave  sequence  thus  contained 
256  data  values.  This  data  digitization  scheme  was  required  In 
order  to  Insure  that  the  data  could  be  FFT  filtered  as  previously 
described  using  an  FFT  algorithm  of  base  2.  The  FFT  algorithm  was 
taken  from  Bloomfield  (1976).  Only  seven  of  the  eight  waves  were 
actually  used  In  the  final  data  analyses. 

The  sonic  profiler  used  to  record  the  Instantaneous  water  ele¬ 
vation  occasionally  suffers  from  data  "drop-outs"  or  "over¬ 
ranges."  These  data  are  Identified  as  "bad-points”  and  are  removed 
numerically  by  linearly  Interpolating  between  the  good  data 
values.  On  rare  occasions,  these  Intervals  of  bad-points  are  ex¬ 
cessively  long  and  result  In  forming  a  fictitious  wave  (see  Wave  #4 
In  Run  #5  in  Fig.  4-5).  These  data  have  been  Ignored  In  the  data 
analyses  because  there  were  only  two  "bad  point"  waves  In  the  total 
of  182  waves  recorded.  These  256  data  points  per  wave  were  reduced 
to  33  data  values  (32  time  Increments)  for  each  Individual  periodic 
wave  In  the  final  data  analyses  for  the  force  coefficients.  The 
total  number  of  discrete  data  values  digitized  for  the  eight  waves 
In  the  41  periodic  runs  that  are  Identified  In  Table  3-13  was 


